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NOMENCLATURE 


mass of rotor i (i = 1,2,3,4) 


pe 
Wy = weight of rotor i (i = 1,2,3,4) 

Jj = mass moment of inertia of rotor i (i = 1,2,3,4) 

xj = displacement of rotor i along x direction(i = 1,2,3,4) 
yj = displacement of rotor i along y direction(i = 1,2/3,4) 
@} = angular displacement of unit i w.r.t. z axis 


(i = 2,2,3,4) 
a, = displacement of tooth spring 
a5 = gear whirl 
phase angle of displacement of unit i (i = 1,2,3,4) 
ej = position of unbalanced mass of unit i from the’ center 
of mass (i = 1,2,3,4) 

= rotating frequency’ of shaft 

= natural frequency of whirl along x direction 
(simple supports with load in the middle) 

= natural frequency of whirl along y direction 

= natural frequency of whirl 

‘0 Nominal natural frequency of gear mesh mode 

‘= uncoupled natural frequency 

= excitation frequency of torque 
$*= period of rotation of shaft 1 (21/e 

period of gear mesh mode (27/4 
period of torque excitation (279u,..) 
period of whirl along y direction f2z/e,, 

period of whirl (27/u,, 4 

= stiffness coefficient of journal bearings of shaft j 
along p direction due to displacement along q 
direction (p = x,y: a= x,y? 3 = 1,2) 

Cpqj = damping coefficient of journal bearings of shaft 5 
along p direction due to displacement along q 
@irection (p = x,y? q= x,y? 5 = 1,2) 

torsional stiffness of shaft 3 (j = 1,2) 

friction coefficient of shaft 3 (4 = 1,2) 

area moment of inertia of shaft j w.r.t. x axis 

Gj = 2,2) 
I,j = area moment of inertia of shaft j w.r.t. z axis 
G = 1,2) 
Ajj = influence coefficient (i = 1,2,3,4; 5 = 1,2,3,4) 
Bi} = coefficient related to bearing performance 
(i = 2,2,3,47 5 = 1,2,3,4) 
= mesh stiffness 
XB = nominal mesh stiffness 
c& = mesh damping 


1) 


vy = mesh frequency 
"= pressure angle 
radius of rotor i (i = 1,2,3,4) 
diameter of shaft j (j ='1,2) 
torque on rotor i (i = 1,2/3,4) 
gear ratio 
width of tooth 
distance of phase plane instantaneous center fron 
origin 
instantaneous radius of arc on phase plane 
instantaneous length of chord on phase plane 
instantaneous angle of arc on phase plane 
Young's modulus 
G = shear modulus 
1, = length of shaft 1 between overhang mass and bearing 
close to it 
L, = length of shaft 1 between gear and bearing 
length of shaft 2 between overhang mass and bearing 
close to it 
Ig, Ug = length of shaft 2 between gear and bearing 
1S t&tal length of shaft 1 (=L, + ly +L) 
¢ = clearance of oil film bearing 
u = relative oil viscosity 
u(t) = oil viscosity at temperature T 
2 = number of teeth 
accr. cr. = predetermined accuracy criterion for 
phase plane calculation 
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Many high speed mechanical systems incorporate gearing 
for speed reduction. This study investigates the stability of 
multi-rotor geared systems supported on oil film bearings 
taking into consideration the coupling between torsional and 
lateral dynamics. 

The emphasis of the study is on the analysis of the 
interaction between the combined torsional and whirl 
instabilities. 

The result of simulation in this study checked the 
previously published empirical approach for predicting 
nonharmonic whirl of geared systems which results from the 


inertial coupling at the gear mesh. 


xii 


The feasibility of inducing a lateral and the torsional 
instability to neutralize an anticipated unstable condition 
is investigated. The possibility of suppressing the 
instabilities by controlling the parameters of the oil film 
bearings is also considered. 

It is interesting to note that the coupling between the 
torsional and lateral oscillations can induce torsional 
instability due to unbalance excitation. Similarly, a 
torsional disturbance can induce whirl. 

A procedure is given to guide the designer in selecting 
the system parameters which insure its safe operation in 
critical situations where instabilities or resonances are 
expected. 

It is shown that torsional instability can be avoided by 
inducing a tuned whirl using an unbalance excitation in a 


system supported on either ball bearings or oil film bearings. 


CHAPTER 1 
INTRODUCTION 


Gears and Geared systems 


Gears, which have been in use for over three thousand 
years, are one of the oldest and most widely used machine 
elements. They proved over the years to be the most important 
means of transmission because of their ability to transform 
wide ranges of power and speed at a high efficiency. 

Historically, gear design has been more art than science, 
but this has changed rapidly in recent years to become an 
important scientific discipline. In their earlier use, the 
main requirement in the design of a geared system was the 
ability to carry the necessary load at the desired speed 
ratio. The load carrying capacity determines the size of the 
tooth, while the speed ratio is satisfied by the relative 
number of teeth in the matching pairs. Because of the early 
slow speed application of gears, no attention was given to the 
instantaneous uniformity of the speed. 

The increase in the speed of operation, however, requires 


consideration of numerous problems such as profile geonetry, 


2 
tooth contact condition, fluctuation of mesh stiffness, 
surface roughness, and oil film thickness. 

In many cases, gear design is the province of a 
specialist. Nonetheless, the design or specification of a gear 


train is only part of overall gear system design. 


Gear System Dynamics 


Gear system dynamics is a problem which has received 
considerable attention due to its practical importance. 

There are numerous sources of excitation in geared 
systems which include errors in manufacturing and assembly as 
well as the time-varying nature of mesh stiffness which is 
inherent to the transmission of motion between the mating 
teeth. These disturbances manifest themselves in torsional as 
well as lateral oscillation which can produce resonances and 
instabilities even under conditions of constant applied loads 
and speeds. 

Instability can be viewed as the result of negative 
damping creating a build-up in the system oscillation. The 
rate of build-up is generally exponential in nature, and 
consequently, very high amplitudes of oscillation can be 
reached in relatively short terms. In the case of geared 
systems, instabilities can occur in the form of rotor whirl 


(lateral oscillation) at certain rotational speeds as well as 
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torsional oscillations at particular mesh and excitation 


frequencies. 


Torsior ics of Gea: 


Torsional vibration of gear systems has been extensively 
treated in the literature [1-21]. It is commonly accepted that 
there are several sources of torsional vibration in geared 
systems. These excitations can be classified as internal 
disturbances which are inherent to the system due to the 
nature of the gear mesh, and external disturbances resulting 
from fluctuations in the input and output loads or speeds. 
Excitations resulting from tooth engagement include the effect 
of tooth error, momentary separation of the teeth, and 
variation of mesh stiffness. 

Houser and Seireg [1] evaluated experimentally the effect 
of gear system parameters on tooth loads at speeds away from 
system resonances. In order to observe the dynamic phenomena 
in gear teeth beyond statistical variations, large pitch 
errors were introduced to insure that the dynamic data were 
sufficiently magnified. Large variations in face width were 
also introduced to simulate significant changes in contact 
patterns. The study shows the result of the test dealing with 
face width variation in both a spur gear and a helical gear 
revealed no noticeable speed effect on the strains developed 


in the gear teeth. The tests involving large manufactured 
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pitch errors showed significant dynamic strain effects which 
appear to be approximately linear with speed and independent 
of load for both spur and helical gears. 

An important source of dynamic excitation in gear systems 
is the variation in mesh stiffness due to the changes in 
length of contact. This change, in the case of spur gears, is 
inherent in the gear operation due to the transmission from 
double to single tooth contact. other sources which contribute 
the variation in mesh stiffness are changes in load 
distribution and change in contact patterns. These may result 
from bearing clearances, shaft and housing deflections, or 
manufacturing errors. The variation of tooth stiffness of a 
pair of matching gears gives rise to a periodic excitation at 
the gear mesh which results in dynamic response of the gear 
system. The oscillatory tooth forces may attain large 
magnitude to cause failure of teeth. Torsional resonance in 
gear systems can occur at frequencies, other than the nominal 
natural frequency, which are called subharmonic and 
superharmonic resonances [2, 3]. 

There have been extensive studies for the prediction of 
regions of system resonance and instability for geared systems 
with variable mesh stiffness [2-5]. 

The accurate evaluation of the actual pattern of 
stiffness variation in a gear system is a difficult problem. 
The compliance of the mesh is influenced by the position of 


the point of contact on the tooth profile, variation in the 


5 
length of contact, shaft misalignment, load distribution, and 
deformation of the wheel body. 

Bollinger and Harker [4] considered the effect of 
harmonically varying mesh stiffness on the system stability 
and calculated the necessary amount of damping for stable 
operation. 

Another important aspect of the torsional dynamics of 
geared systems is the resonant behavior of the pinion-gear 
system under steady state operating conditions. 

Benton and Seireg [2] presented a computer simulation 
procedure for predicting the steady state response, 
resonances, and instabilities of pinion-gear system subjected 
to sinusoidal torsional excitation. Their study isolated the 
gear pair from the rest of the system and treated it as 
uncoupled from the other connected elements of the total drive 
based on the assumption that mesh of gear pair has high 
stiffness and low effective inertias. Sinusoidal and square 
wave variations of the mesh stiffness were considered to 
investigate the dynamics of the system under different speeds 
when subjected to uniform and cyclic loading conditions. A 
response spectrum of a simple model of a pinion-gear pair 
subjected to harmonic excitation was presented. The results 
which were correlated with experimental data show the 
influence of the harmonic excitation on the overall response 
of the system. The different models of stiffness variation 


used showed considerable influence on the dynamic behavior of 
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pinion-gear systems. Variations of mesh stiffness affects both 
the frequencies and the pattern of the geared system 
resonances and instabilities. 

Benton and Seireg [3] extended their research on 
torsional vibration of pinion-gear system to investigate the 
parameters causing high dynamic magnifications and 
instabilities in multiple reduction geared systems. These 
factors include the effect of external excitations, system 
inertia, variation in mesh stiffness, contact ratio, and 
damping in the mesh. The study presented a procedure which 
can be applied to multiple reduction gear drives to aid in 
the analysis and design of such systems when the exact nature 
of stiffness variation is not precisely known. The influence 
of the system inertia, variation in mesh stiffness, contact 
ratio, and dissipation in each mesh on the gear system 
performance is also discussed. 

Tooth backlash can also be a contributing factor in the 
vibrations of geared systems. Azar and Crossley [5] examined 
these effects both experimentally and analytically in an 
attempt to predict the vibrations observed in the driven shaft 
of a lightly loaded gear system. They reported the presence 
of large vibrations which are influenced by the tooth backlash 
in the driven shaft when the mesh frequency is one half the 
torsional natural frequency of the driven shaft system. A 
digital simulation model, in the study, is developed to 


represent a lightly loaded geared torsional system consisting 
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of a drive unit, gear pair, and load connected by flexible 
shafts. 

Problem of tooth separation and impact in a geared system 
has also been subjected to considerable attention (5-10). 
Innoue and Seireg [6] presented a digital simulation of the 
vibration and the impact in a dual drive branched gear system 
during the engagement and disengagement of the prime mover. 
The study illustrated the existence of an optimal backlash 
condition for the idle pinion which minimizes the dynamic 
tooth load. Seireg, Shah, and Khazekhan [7] investigated the 
effect of gear vibration through the backlash on the dynamic 
tooth load experimentally. 

other factors which place a penalty on the load carrying 
capacity of gears are the manufacturing errors. The types of 
errors and their effects on the system dynamics can be found 
in previous works [11-13]. 

‘There are many practical situations where resonances and 
instabilities in pinion-gear systems are difficult to predict 
in the design stage due to the unreliability of estimating the 
mesh stiffness and damping parameters. 

In an attempt to insure the stability and safe 
performance of pinion-gear systems, Benton and Seireg [14] 
investigated the feasibility of the design of vibration 
absorber for gear drives which can provide an appropriate 
solution for such conditions. They developed the dynamic 


equations of motion for the absorber system to determine the 
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pertinent variables necessary for general description of the 
absorber system. The conditions for uncoupling these equations 
are determined, and the design parameters to be considered in 
the optimization are specified. A search is undertaken to 
evaluate the optimal parameters for different gear system 
conditions and the results are given as design charts. The 
study shows that an optimal dynamic absorber can be designed 
which is capable of providing safe operation of pinion-gear 
system at all speeds including unexpected resonance and 
instability conditions. 


tera amics of Gear Systems 


A considerable number of studies of gear system whirl 
have also been reported. A review of the problem and 
developments in the theory of whirl can be found in previous 
works [22-44]. Some of the recent contributions to the 
analysis of whirl are the effects of gear constraint on the 
whirl of simple supported pinion shafts [22], the effects of 
the change in the cross section [23-25], and the effects of 
shaft support (26, 27). 

It was reported that the effect of the gear restraint on 
the whirl of the pinion shaft is a reduction of the whirl 
amplitude, a change in the resonant frequency, and the 
appearance of subharmonic as well as superharmonic resonances 


(22). 
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A simulation of the whirl interaction at the gear mesh 
is investigated, in which the pinion-gear system is modeled 
as a coupled two rotor problem [28]. The gear rotor, in the 
study, is assumed to be supported on isoviscous fluid bearings 
and the transmitted gear load is taken into consideration. 
Each rotor is simulated as a two degree of freedom system 
coupled at the point of contact during the mesh. The result 
of the study shows that the gear constraint with zero 
transmitted load causes an increase in the threshold speed. 
Instability in this case occurs, as expected, in both the 
pinion and the gear rotors. The effect of the transmitted load 
is also shown to increase the threshold speed with one of the 
two rotors reaching instability first and therefore causing 
the system instability. 

A rotor whirl simulation procedure [29] is described for 
investigating the effect of the unbalance magnitude on the 
whirl of rotors. In the study, a rigid rotor on a simplified 
isoviscous film is considered for illustration in an attempt 
to present a general simulation procedure which can be used 
for the analysis of this class of problems. 

The designer is generally interested in predicting the 
expected peak amplitude of the whirl orbit at a given 
operating speed resulting from a particular amount of 
unbalance as well as the stability of the rotor operation in 
the design stage. Some of the contributions to the analysis 


of this problem are made by Lund et al [30-38]. 
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The effect of fluid film bearing on the whirl has also 
been widely studied. Seireg and Dandage [39] described an 
empirical design procedure for predicting the performance of 
fluid film bearing based on experimental thermohydrodynamic 
consideration. The study utilized developed relationships to 
construct design nomograms for evaluating a modified 
Sommerfeld number which can replace the isoviscous Sommerfeld 
number in evaluating the performance characteristics of the 
bearing for any given oil and inlet temperature. The modified 
Sommerfeld number can then be utilized in the standard 
formulas to calculate eccentricity ratios, oil flow, 
frictional loss, and temperature rise, as well as stiffness 
and damping coefficients for bearings. 

It is generally accepted that the dynamic properties of 
the fluid film can be expressed in terms of four stiffness 
and four damping coefficients (39, 45, 46]. The investigation 
of rotor stability involves the computation of the stiffness 


and the damping coefficients. 


Coupling of the Lateral and Torsional pynamics 


It has been common practice for the design procedure of 
@ geared system to analyze the torsional dynamics and the 
lateral dynamics separately. In this way, the effect of gear 
coupling is only partially accounted for, and the coupling 


between torsional and lateral vibrations is neglected. As a 


un 
result, the procedure fails to account for several important 
aspects such as the shift in critical speeds and change in 
mode shapes caused by the coupling, change in the instability 
threshold, the excitation of all rotors in the system caused 
by the unbalance in any of the rotors, and the damping in the 
torsional modes contributed by the lateral motion of the 
journals in their bearings. 

Lund [47] used the impedance matching method [16] in an 
effort to calculate the coupled torsional-lateral vibrations 
in a geared system assuming a constant mesh stiffness. The 
study considered the damped forced vibrations of the system 
caused by mesh errors and unbalance, whose complex eigenvalues 
define the damped critical speeds and the instability 
threshold of the system. The vibratory modes for the rotors 
supported on fluid film bearings are calculated independently, 
using the Holzer method for torsional vibrations, and the 
Myklestad-Prohl method for lateral vibrations. They are 
subsequently coupled through impedance matching at the gear 
meshes. The resulting equations are solved for the unknown 
mesh contact forces, and the roots of the coefficient matrix 


give the eigenvalues of the system. 


Scope of This Work 


It can be seen from the literature review that most of 


the previous studies dealt with either the lateral or 
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torsional dynamics separately for isolated pinion-gear 
systems. Even in the studies treating the entire geared 
systems, the investigations of dynamics of the system were 
confined to specific phenomena. Effects of coupling of 
torsional-lateral dynamics were treated by Lund [43]. He 
considered the coupling with the assumption of constant 
stiffness for the mesh in order to simplify the determination 
of the coupled frequencies of the system and the threshold of 
instabilities. 

Since the variation of stiffness during the mesh cycle 
is an inherent characteristic of gear drives, this study is 
therefore initiated to investigate the combined effects of 
coupling the torsional vibration with the consideration of 
time-varying stiffness and lateral dynamics of a multi-mass 
geared system. The drive shafts are supported on fluid film 
bearings and the rotors are subjected to unbalance and torque 
excitations which can be expected in many practical 
applications. An important objective is to develop rules which 
can be followed in the design stage in order to insure the 
safe operation at all operating conditions including 
resonances and instabilities. Emphasis is placed on the 
interaction at resonant and unstable regions. The goal is to 
develop a procedure for selecting the system parameters which 
insure the stability of the system by taking advantage of the 


coupling effects. 
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Since both the lateral and torsional oscillations of a 
geared system interact at the points of the tooth contact, it 
is conceptually possible that a system can be designed where 
the oscillations at the contact point can be made to 
counteract and vanish by the appropriate selection of the 
design parameters and the phase relation between the 
excitations. 

A simple supported geared system with two overhang rotors 
is considered in this investigation. Each rotor has three 
degrees of freedom. As a result, the system has twelve degrees 
of freedom. In unstable situations, the response of the system 
usually exhibits one dominant frequency of oscillation at 
which build-up occurs. 

‘The possibility of suppressing an instability by inducing 
another instability of a same type is investigated. The basic 
concept is extended to see if it is also possible to suppress 
an instability by inducing another "tuned" instability of a 
different type. Since the system under consideration has time- 
varying coefficients and is supported on bearings with non- 
linear characteristics, exact analytical stability criteria 
are not practical. A procedure based on numerical solution of 
the equations of motion is developed with the objective of 
guiding the selection of the system parameters which insure 
its stable performance under the expected operating 


conditions. Because of the numerical nature of the solution, 
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the developed approach can be readily extended to the design 


of more complex gear systems. 


CHAPTER 2 
SYSTEM SIMULATION 


hysical Model 


The systems dealt with in this work are typical geared 
systems. They represent a single reduction gear unit with and 
without overhang rotors. 

The systems, which are illustrated in Fig. 2.1 and Fig. 
2.2, respectively, are composed of a gear pair with the high 
and low speed shafts simply supported on fluid film bearings. 
The shafts are assumed to be massless. Each of the gear wheels 
and the overhang masses is assumed to have three degrees of 
freedom representing the x and y movements of the center of 
mass and the angular oscillation, © about the center axis. 
Since the analysis is primarily directed toward high speed 
systems, the outer diameters of the wheels are relatively 
small and consequently the gyroscopic movement of the gears 
can be neglected. The system parameters are listed in the 


Nomenclature. 
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Gear (rotor +) 


2 ie wy 


cine 


Fig. 2.1 Two Rotor Geared system 
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drive (tor 2) gear (rotor ) 
ye : “$s y 


pinion (rotor 3) load (rotor +) 


Fig. 2.2 Four Rotor Geared system 


ae 


uy 


1s 


syst ons 


The equations of motion can be derived in the same way 
for both systems. Therefore, the equations of motion for the 
four rotor geared system are obtained, and the two rotor 
system is considered to be a special case of the general 
formulation where the overhang masses are zero. Since the 
general system has twelve degrees of freedom, a 12 x 12 matrix 
(Fig. 2.3) consisting of twelve equations of twelve variables 
can be constructed as follows. 

The coefficients in the matrix are obtained by assuming 
a small movement in each degree of freedom and evaluating the 
forces induced by this movement on the entire system. For 
example, for element 1, a small positive displacement x, is 
imposed on the system, and the resultant force on rotor 1 
along the x direction is obtained to fill that element. The 
resultant force is evaluated from the reactions at the mesh, 
the shafts, and the bearings as a response to the movement. 
For the rest of the elements in the first column, a similar 
procedure is applied to obtain the forces on all four rotors. 
Also, the remaining elements in the matrix are determined by 


evaluating the forces due to the displacements x2, x3, -+ 


4. 
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Fig. 2.3 A 12x12 Matrix of Coefficients for the 
Equations of Motion 
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eat tl sh 
Fig. 2.4 illustrates the geometry of a pair of gears [48] 
in mesh. The force induced on rotor 1 due to the displacement 
X, can be expressed as 
ax 
0S [ky (-cosd) x) + cy (-cosp) —* } 2.2 
at 
where 


k, , the time varying stiffness of the 
mesh 


= Ky (1 - A cos (wy t)) es 


(kg: the nominal mesh stiffness) 


¢, , the damping coefficient at the mesh 


My ot 
2€ o —+—}- 2.3 
My +My 


, , the mesh frequency 
= (rotating speed of shaft) 
(number of teeth of gear) 2.4 
A, constant representing the degree 


of fluctuation of stiffness 


= .5 for spur gears 2.5 


€ , damping ratio 


(assumed to be 0) 2.6 


@, pressure angle 


(assumed to be 20 degree) a7 
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Force at the Bearings 

A diagrammatic illustration of the effect of the oil 
film bearing support, and the motion in the y-z plane are 
shown in Fig. 2.5. 

The stiffness and the damping coefficients due to the 
bearing support is calculated using the subroutine Brng 
(Appendix) which is based on the work by Seireg and Dandage 


[39]. The force on rotor 1 can be written as 


Force due to Bending 
As illustrated in Fig. 2.6, the force due to the bending 
of the shaft [49] as a result of displacement x, can be 


expressed as 


6 ET, (Ly + by) 
2 


Ly [Ly? + Lg? - (Ly + L5)?) 


ka = ky [1 — A cos(w, t)] 


Kp, : My Mg 


c=2 cH) [ky s 


1 


Wye My 


Fig. 2.4 Stiffness and Damping from Tooth Contact 
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(Ly tLe +L) 


(lL. +L.) 


Fig. 2.5 Shaft Supported on Oil Film Bearings 


23 


24 


a a 
7] 7 z 
F LeLolle't Le'- et Lo] 
= sri 


6ETetLs) 


Fig. 2.6 Shaft Subjected to Bending 
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Equations of Motion 
The 144 elements of the matrix (Fig. 2.3) can then be 
used to formulate the following equations of motion for the 


general system. 


2, 
%1 (C0876 Ky + Kye Bay + 


A, 


11. 
ax. 


(©0876 Gq + Cyoey Byy) 


at 


~ ¥, (cose sing km + Kyyq By) 
ay. 
at 


(cose Sind Cy + Cyyq By) 


ao. 


- ©; r, cose ky - ry, cosd cy 


at 
1 ax, 
“Pitty Osa Bag St rae Be 
Aya at 
ay, 
+ ¥2 My Bia + F, Sxya Baa 
at 
ax. 
+ x, cos"g km + — cos”¢ cy 
at 
ay. 
+ y3 cosp sing ky + cosp sing Cy 
at 
ao. 
~ ©; ry cos ky - —* rq cos¢ «, 
Ee) Pa ae, m 


+My ey ¥,? cos(a, t + %) 


2 cos(w, t + $3) 2.10 
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x (CosP Sind ky + Ky By) 


ax, 
Fe (0088 SiN Sy + Cyr Baa) 


1 

sn? 

Ya (i Kg + Kya Big t+ 2) 
fet Man Pat 5 

ey, 

at 


tn’, 
(sin?@ co, + ey, Bay) 


ae. 


6, ry sing k, - 
zt Laas 


2 Kya Bia + 


1 ay, 
yp ( By +—) +—*> ey, B. 
2 yya Bre Rg dee ara eae 
ax, 
x, Sin@ cosd ky + sing cosp c, 
at 
ay. 
in? a) ad, 
ys sin?o ky + sin?9 «, 
i m 
ao. 
3X3 sing k, - — ry sing cy 
at 
My e, 0? sin(w, t + #4) 
a, 
22 2 


sin(w, t + ®) 


#3 


ax. 

Sp Som ey 
F ayy 

Eelam 


2 de. 
Oren sista halo 


Key + X34, COSd ky + 


ay. 


r, sing i r. 
1 Sind yt Fry 


1 a 
Ls © a a ar 
Aya ai 
hay yg 
+o 
aya Pia Se Saye 


cos? Cy 


sing oy 


x4", 


ax. 


Ty COS# Cy 


sing cy 


x. 


Spy 8B: 
a aa 
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ax. 


Kya Bap #2 
1 Kya at 


a ay, 


+ ya Oya Ba #) t 


A, at 


12 


ax, 


Per ae 
2 %pa Ba 7 3 


= Yo Oy, Baa + 
‘22 


+My e2 W,” sin(w t + #2) 


Aa 


Syyr Pia 


+ > my ey 0,7 sin(w, t + %) 


Aye 


x, cos’g k, + 


at 


ay. 


+ y, cosp sing k, + 


ae. 


+ 1 Fy C086 hy + 


2, 
X3 (€087b ky + Kyyo By3 + 


A 


cos sing cy 


Fy COSd Cy 
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ax, 


= (00876 cq + Cyyp By) 


~ ¥3 (C08 Sind ky + Kyyp By3) 


ey5 P 
- —=— (cos sing cy + c, 
at 


‘xy2 B33) 


cosp ky + 3 COsp c, 


+ 4 aca Bag + 


¥q Kgpa Bag 
4 Kaye Baa + ob 


Sxy2 Baa 


+My e3 0,7 cos(w, t + #3) 
A 
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+S my eg 62? cos(a, t + 04) 
Asa 


ax, 
cosp sing ky + cose sing cy 
ae 
ay 
+, sin’s k, + —* sin’ o, 


at 


+ ©, r, sing ky + 


1 Sind Cy 


at 


= X3 (COS® Sind ky + Kyys B33) 


ax. 


Fe (SOSH SiMe Cy + Cyyy Byq) 
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1 
~ ¥3 (sin?g k, + Kyy2 By3 + =) 


+ 0, ry sing k, + 


ax, 
* 4 Kyxa Bag + 2 Syn Bag 
at 
a ay, 
Peta: Oyya Byars Be 
agg at 


sin(w, t + #3) 


4 M2? sin(u, t + %4) 


ax. 
~ 1 Fy C088 hy - > ry Coss Cy 

ayy ri 
~¥1 Fy sind hy - ry sind cy 

a0, 

Pr RaEa fa Maisie = Ea iG 

ax, 
+ 5 Ty CORD ky + TF ry C088 Cy 

ay, 
+ ¥5 Fy sin ky + > ry sino cy 
+8, 
+8 Dy 2.18 


¥3 Kyxa 


%4 OSoa Bag + 


¥3 (Kyy2 


9 Kyxa Bag + 


By, + 2 ». 
+o, 
24 + Syne Baa 
1 ax, 
ri 
= ce 8 
wea Pas 
Mae ae 
B, Bs 
¢, 
a4 5 Sxya Bae 


cos(w, t + 65) 


ax. 


it 


Byq + 


4 Myxa Bag > 
oor 

$y beats oe 

4 ya Baa +) ~ A eye Bag 

My e4 7 sin(w, t + 04) 


2° sin(w, t + &3) 


a ae, 
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(es eae 2 


at 


Bees 


a1 


2.21 


32 


tical Program stru 


The equations obtained in the previous section are 
Programmed with the structure detailed in Fig. 2.9 and Fig. 
2-10. The main technique used to solve twelve coupled, time 
varying differential equations is the phase plane method which 
is a versatile tool for the analysis of this type of systen. 

Twelve phase planes are monitored simultaneously to 
analyze the response of the system shown in Fig. 2.2. The 
primary advantage of the technique over other numerical 
technique is that integration is based on the incremental 
angie in the phase plane rather than time. This allows very 
efficient integration of the equation of motion for systens 
with a wide band response spectrum. 

The phase plane for undamped free vibration of a single 
degree of freedom system is shown in Fig. 2.7. The motion has 
an amplitude A, natural frequency ,, and the phase angle y. 
The displacement and the velocity can be reduced to a single 


equation representing a circle. 


2.22 


This method can be applied to complicated systems by 
rewriting the equations of motion as in Fig. 2.8. The 
integration procedure is shown in Fig. 2.9, and Pig. 2.10. 
Details can be found in the previous works by Seireg et al 


(29, 50). 
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Fig. 2.7 Phase Plane for Undamped Free Vibration 
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Fig. 2.8 Phase Plane § Method 
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Given: geometry, oil temp., r-p.m., 
material constant, torque, 
avg. mesh stiffness, G.R., 
initial conditions 


if 
call Brng to calculate 
stiffness and damping coefficient 
I 
I<o 


>—————4 


ue 


Integration step for one variable 
with incremental angle related to 
its eigenvalue 


Repeat the steps shown in Fig. 2.10 
for the second variable 


Repeat the steps shown in Fig. 2.10 
for the last variable 


ye! < Tag 


no 


Fig. 2.9 Structure of Program to Solve Equations of Motion 
Using Phase Plane 6 Method 
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A 


Start the integration step 
by phase plane § method 


I+ 0 
a] 
Calculate 6 
ax 4 
a tpt (x + sj-0 
at’ 


Calculate the instantaneous radius 
of arc and the length of chord 
ax. 


= LO + 53)? + ¢ 


Pat 


3 sin( 


) 


2 
I 


Call Psis(to decide 9;) 


Call Phase(to define the next position 
on the phase plane) 

Pp dt 

xy £1; sin(a; + ra) 


Mia 
ax dx 


(Vier = (yj 
aes tt att 


Fig. 2.10 Details of Integration step in Fig. 2.9 
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Iteration 


(use Hi = ED py, 


to calculate 6;, 1;, 


e am 
and (—) 
ier and 


it) 


and obtain 


Calculate Ax 
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Bee Ot. Ra)? + 0 ga > 


Fas pacer: = 
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G41 ita no 


tig = ty + at 
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Fig. 2.10 continued 
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Significant et rs ste! 


Since all the unstable situations are related to the 
natural frequencies of the system, it is helpful to calculate 
the uncoupled natural frequencies of individual modes for the 
isolated system to have a better insight into the coupled 
system. Geared systems, such as those investigated in this 
study, have many uncoupled modes as follows. 

sion of Sha: 

The natural frequency of the torsional vibration due to 

the torsional stiffness of shaft between rotor 1 and rotor 2 


(or rotor 3 and rotor 4) is 


ex 
Ne 
“5m Fy? +m) 
coar_systen 


Assuming that the motion is confined along the line of 
action, the natural frequency of gear-pinion mesh (Fig. 2.4) 


can be written as 
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Lateral Vibration of Simple Supported Mass 
due to Shaft Bending 


For this mode, assuming a centrally located mass, the 


natural frequency can be calculated as follows: 


where 


te: ue 


The natural frequency in this case can be expressed as 


k, 
bo 2.27 
a 
where 
6 EI 
pee ; 2 2.28 
2 Ly? (Ly + Ly + Ly) 64 


The system, investigated in this study, is assumed to be 


symmetrical for simplification. 


CHAPTER 3 
SYSTEM RESPONSE AND STABILITY CONSIDERATION 


The system investigated in the reported study is assumed 
to be symmetrical in order to simplify the computation. The 
following assumptions are made: M, = My, Ly = Ly = Ly = Lg: 
Several cases are investigated to study the effect of the 
design parameters on the system stability and the possibility 
of introducing additional instability to neutralize the 
adverse effect of an expected instability. 

In geared systems, instabilities can occur in the forn 
of the lateral as well as torsional oscillations. The whirling 
and the torsional motion interact at the points of tooth 
contact. Accordingly, the amplitude build-ups can possibly be 
made to counteract and vanish. 

Unstable situations are selected from the regions which 
represent practical operating speeds, and each case is 
stabilized by an additional excitation and by proper selection 
of the design and operating parameters including the bearing 
geometry, oils, geometry of the gear, and the phase shift of 
the induced disturbances. 

Although the time history of the deflection of tooth, 


de, and gear whirl, dj, are selected to represent the 
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torsional and the lateral dynamics of the system respectively, 
all the twelve main variables(x,, x2, +», Yy, - 


calculated simultaneously and are examined when needed. The 


» @4) are 


values of dy and a, are computed as 


4 


The following are the main system parameters used to 


investigate the system used in the reported study. 


@ = 20 degree 

0, = 21 N/ 60 rad/sec 

01 (GR) rad/sec 

Om = 2 (# Of teeth of pinion) rad/sec 
Dy = 2.5" 

Dy = 2.5" 

d,,(bearing diameter) = D + 2.¢ inch 
1, (bearing length) = 2.5 

Lg = bg = 4.325" 

+169 slug 


Wy (1000) lb-in 
ry 1, /GR 

Kg = 1-2 (10)6 ibsin 
E = 30 (10°) 1b/in? 
Poisson's ratio = .292 
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The Selection of Bearing Parameters 


The parameters of the oil film bearing for different 
situations are calculated as follows. In the case of stable 
system conditions, any bearing parameters which do not create 
negative damping are acceptable. However, in the case of 
system instability, the bearing parameters can be selected to 
supply the necessary amount of positive damping to the systen. 
The system parameters are used to calculate the bearing 
parameters for the two rotor system. In the reported study, 
SAE 20 oil is selected as lubricant whose relative viscosity, 
u, is 1.36x107® reyns. Oil inlet temperature is assumed to be 
110° F, and iteration method [39] is used to calculate the 
temperature increase, and accordingly viscosities at those 
operating temperatures which are calculated from the 


relationship 


M(T) = & exp(1271.5/(T + 95)) 


along with the Sommerfeld numbers. Then, the Sommerfeld 
numbers are used to calculate the stiffness and the damping 
coefficients for given systems in different operating 
conditions. 

Table 3.1 shows radial clearance of the fluid film 
bearings for each operating speed which is selected in an 


iterative way to reach the maximum damping for given 


Table 3.1 Selection of Oil Film Bearings 


Speed (rpm) Radial Clearance (inch: 
4200 -002 
23875 +004 
35812 -005 
47750 -006 
95500 -006 
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44 
conditions. The subprogram to calculate the stiffness and the 


damping coefficients is listed in the Appendix. 


or System w. ts 


A two rotor system with simply supported shafts is 
considered in this section. Three cases which represent 
stable, resonant, and unstable operations [2] are considered 
here to investigate the effect of using additional torque 
excitations and oil film bearings to suppress the amplitude 
of the oscillation due to a fluctuating torque input. 

The system used in this section is shown in Fig. 2.1. It 
has six degrees of freedom and is originally supported on ball 
bearings. The system parameters are as mentioned before. 

table Conditions 

In this section a case is selected where no build-up 
occurs in the system response. If the amplitude of the 
response is acceptable, no additional treatment is necessary. 
However, if the response amplitude is considered to be 
excessive, it can be reduced by an additional torque 
excitation and an appropriate design of the oil film bearing. 

0 To! tation on 

The stable response of a two degree of freedom torsional 
system due to an oscillating torque input is calculated for 
the following conditions for the same system shown in Fig. 


22s 
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Ty = Wy (1000) 1b-in 
gy = 3750 rad/sec 
Uno = 5012 rad/sec 


3750 rad/sec 


Fig. 3.1 shows the deflection at the tooth contact. a 
beat can be seen in the response, but no build-up occurs in 
the motion. The beat is attributed to the difference between 
natural frequency and the excitation frequency. The frequency 
of oscillation is the same as the excitation frequency, and 
the beat frequency is approximately 1250 rad/sec which is 
equal to the difference between the natural frequency and the 


excitation frequency. 


of an onal Torqui ut_o1 

An additional torque input which can give the same 
response amplitude is now imposed on the other gear. As a 
result, the amplitude of the oscillation becomes much smaller 
as shown in Fig. 3.2. Elimination of the large amplitude is 
essentially complete. Practically, there is no deflection of 
tooth spring. (The peak amplitude is reduced from 0.0004" to 


0. 0000006" 


In this case the amplitude of the second torque 
excitation is the same as the first torque excitation and the 


Phase angle between the two is 180 degrees. 
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Oey = 3750 rad/sec 


Case 3: The Effect of the Oil Film Bearings 


Fluid film bearings with appropriate parameters (oil 
viscosity, clearance) is used to replace the original ball 
bearings. The damping in the bearing is utilized to reduce 
the amplitude of the oscillation of the system with the two 
torque inputs. The bearing length, in this study, is fixed 
with the ratio of 1,/d, = 1.0. For this case, the following 


parameter values are selected. 


© = 0.005" 


Be 1.36/(10)® reyns 


As shown in Fig. 3.3, the amplitude of the tooth spring 
response is slightly reduced by the use of the oil film 
bearings. It can also be seen that the beat has disappeared 
as a result of the damping in the bearings. 

It is interesting to note that as a result of using the 
oil film bearing, the system experiences a larger amount of 
whirling of the gear which was negligible when ball bearings 
were used. The amplitude of the whirl, however, is 
approximately 0.0013 % of the radial clearance and can 


therefore be tolerated (Fig. 3.4). 


Fig. 3.1 


00040 pa YY 
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Oscillatory Deflection of the Teeth at the Mesh due 
to a Torsional Excitation at Rotor 1 (Case 1) 
Tex = +0017, The = .00125 
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Fig. 3.2 Oscillatory Deflection of the Teeth at the 
Mesh due to Two Torques (Another Identical 
Torque with 180 Degree Phase Shift Imposed 
on the Other Rotor (Case 2) 
Tex = -0017, 7, = 00125 
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-0.0000010-+ 


-0.0000020 1 1 


49 


0.0000020— 


0.00000104 


ac AU RV er 


9000 0020 


Fig. 3.3 Tooth Deflection at the Mesh due to Two 
Torques When the Shafts are Supported by 
Appropriate Oil Film Bearings (Case 3) 
Tex = -0017, Thy = 00125 
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Fig. 3.4 Gear Whirl for Case 3 
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Resonance condition: 
The system experiences resonance when the exciting 
frequency is the same as the nominal natural frequency of the 


pinion-gear system. 


ster an 

The resonant response of the tooth spring displacement 
in the undamped condition due to an oscillatory torque 
excitation is expected to build up linearly with time. 

Fig. 3.5 shows the response which is obtained for the 


following conditions. 


T, = W, (1000) 1b-in 


© 
0 


= 5000 rad/sec 


© 
" 


5012 rad/sec 


= 5000 rad/sec 


‘The frequency of the response is the same as the nominal 


natural frequency of the pinion-gear system. 
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fase 5: The Effect of the Additional Torque Input 
Another torque excitation at the resonant frequency is 
imposed on the other mass in order to neutralize the effect 
of the resonance resulting from the first excitation. 
In this case the amount of the additional torque is same 
as the original one. A 180 degrees of phase shift of torque 


excitation was found to produce the best results. 


- 7 


5000 rad/sec 


The result which is plotted in Fig. 3.6 shows that the 
amplitude of the oscillation is suppressed drastically and 
the amplitudes of oscillation is approximately 0.1 % of the 


original value. 


se 6 ct_of the 0. 

Bearing length, oil viscosity, and clearance of the 
bearing can be selected in such a way to provide positive 
damping to the system with two torque excitations. For this 


case, the following parameter values are selected. 


0 
nl 


0.006" 


1.36/(10)® reyns 


e 
" 
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Fig. 3.5 Tooth Deflection at the Mesh for a Resonance 
due to a Torque Applied on Rotor 1 (Case 4) 
Tex = -00125, Th, = .00125 
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Fig. 3.6 Tooth Deflection at the Mesh due to Two 


Torque Applied on Two Rotors at Resonant 
Frequency (Case 5) 
Tex = -00125, Tho = .00125 
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Fig. 3.7 Response due to Two Torque Excitations at 
Resonance with Appropriate Oil Film Bearing 
Support (Case 6) 
Tex = -00125, Th. = .00125 
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Fig. 3.8 Gear Whirl for case 6 
Tex = -00125, Th4 = .00125 
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As shown in Fig. 3.7, the response of the system has been 
suppressed more by the positive damping from the selected oil 
film bearing. 

As in the stable case, oil film bearing causes gear whirl 
as shown in Fig. 3.8, which is approximately 0.0013 % of the 
clearance and is therefore of an acceptable magnitude. 

Uns sys’ 

An unstable situation, characterized by exponential 

build-up of the response amplitude, is expected under the 


following condition 


This condition is selected based on the data given in 


Fig. 7(a) of the work by Benton and Seireg [2]. 


onse of the Uns! 

The unstable response is obtained from condition which 
satisfies the relationship expressed in Eq. 3.4 between the 
excitation frequency, the nominal natural frequency of the 


pinion-gear system, and the mesh frequency. 


T, = Wy (1000) 1b-in 


Wey = 2500 rad/sec. 
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Wp = 5012 rad/sec 
ig = 10000 rad/sec 


6, = 2500 rad/sec 


As shown in Fig. 3.9, torsional instability of the system 
is characterized by the exponential build-up pattern at the 
nominal natural frequency of the system which is twice the 
excitation frequency. This type of instability is called 
superharmonic resonance. 


case 


‘The of al ‘orque 

As in the resonance case, an additional torque input with 
identical system parameters and the proper phase shift can 
neutralize the effect of the instability due to an existing 
torque. 

The response at the mesh with two torque inputs is shown 
in Fig. 3.10. The amplitude of the response has been reduced 
approximately to 0.1 % of the amplitude of the response due 
to one torque excitation by imposing the additional torque 


excitation with 180 degrees of phase shift. 


ahaa 


Wg, = 2500 rad/sec 
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fase 9: The Effect of the 0i1 Film Bearings 
As in the example considered in case 6, the resonance 
amplitudes can be further reduced by using fluid film bearings 
to support the shafts. 
The bearing parameters are selected as in previous cases 
to further reduce the torsional oscillations without inducing 
significant whirl. For this case, the following parameter 


values are selected. 


© = 0.004" 


# = 1.367(10)° reyns 


As shown in Fig. 3.11, the amplitude of the response is 
Suppressed and the build-up no longer exists. 

‘The dominant frequency of the response is now the same as 
the excitation frequency, since the superharmonic resonance 
@ue to the negative damping created in the systen is 
Suppressed by the positive damping from the appropriately 
selected oil film bearing. The beat which appears in the 
response is due to the difference between the excitation 
frequency and the natural frequency. As in the previous cases, 
fhe oil film causes lateral gear oscillation as shown in Fig. 
3-12. Again, the amplitude of gear whirl is less than 0.01 ¢ 
of the radial clearance and is accordingly an acceptable 


magnitude. 
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Fig. 3.9 System Response at the Tooth Mesh for a 
Superharmonic Instability (case 7) 
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Fig. 3.10 Response due to Two Torque Excitations at 
Superharmonic Resonance (Case 8) 
Tey = +0025, Th, = .00125 
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Fig. 3.11 Response due to Two Torque Excitations at 
Superharmonic Resonance with Appropriate 
Oil Film Bearing Support (case 9) 
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Fig. 3.12 Gear Whirl for case 9 


Tex = :0025, Th, = .00125 
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Tt should be noted here that in this case the response 

of the system supported on the selected oil film bearings due 

fo each torque input shows the same pattern of vibration. Two 

ferque inputs, therefore, can neutralize the effect of each 
other with the proper amount of phase shift. 

in other words, the large amplitude of oscillation due 

*o a torque input can be avoided by another torque input with 

‘he proper amount of phase shift. The response of the system 

can be further suppressed by the damping provided by 

appropriate selection of the parameters of the oil film 


bearings. 


wi ce Excitat: 


The case of the coupled two rotor system subject to whirl 
is investigated in this section. 

The effects of the gear coupling on the whirl of the 
system are a reduction of the whirl amplitude, a change in 
the resonant frequency, and the appearance of subharmonic as 
well as superharmonic resonances (22). 

The possibility of suppressing the whirl due to an 
unbalance of one rotor by introducing an additional unbalance 
on the other rotor will be investigated for stable, resonant, 


and subharmonic resonant conditions. 
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The system is the same as that considered in section 3.1. 
The shafts are simply supported on ball bearings and the 
system has six degrees of freedom. 
Stable Condition 

4 stable system without build-up of the whirl amplitude 
is investigated here. Even without build-up, the response of 
the system can reach unacceptable magnitudes if large amounts 
ef unbalance exist. Either the amount of unbalance can be 
reduced or another unbalance can be imposed on the other gear 
to suppress the system response to the allowable limit. The 
latter case is discussed in this section in order to deal with 
high speed gear situations where perfect balance is difficult 
to achieve. 
Case 10: Stable Response 

The stable response of the coupled system due to an 
unbalance on one rotor is obtained for the same system 
dilustrated in Fig. 2.1. The following conditions are 


considered. 


+00003/M, inch 


3750 rad/sec 


y = 5029 rad/sec 


The whirl frequency can be calculated as follows. 


2) 4-5 
My = 5 (4? + 9y?)] 3.5 
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where 


3.6 
_ 4B a, ey 
—— a 
(Ly + Ly) 


‘The value of the equivalent natural frequency w, for the 
Geared system is calculated according to the procedure given 
in the previous work by Seireg [20]. The lateral natural 


frequency along y direction, is calculated without 


considering the gear constraint. The frequency 0, is the 
natural frequency calculated in the same manner as oy with the 
exception of adding to the mass of the pinion an equivalent 
gear mass obtained by dividing the moment of inertia of the 
gear by the square of its pitch radius. 

Since the system analyzed in this section is coupled at 
the mesh, an induced torsional motion is to be expected. The 
amplitude, however, of the torsional motion due to the 
unbalance is much less than that of gear whirl. Accordingly, 
the work will be focused on reducing the amplitude of the gear 
whirl. 

‘The time history of the gear whirl is shown in Fig. 3.13. 
#£ can be seen that there is no build-up in this case as would 
be expected. The frequency of the response is the same as the 


rotation frequency of the shaft. 
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case Effect o! ance 

An additional unbalance producing the same amount of 
eccentricity on the other rotor with 0 degree phase shift is 
introduced in the system to suppress the oscillation caused 
by the first unbalance. The result is illustrated in Fig. 
3.14, which shows that the oscillation of the gear is 
essentially eliminated by introducing the additional unbalance 


on the other rotor. 


The parameters of the oil film bearing considered in 


this analysis are as follows: 


© = .005" 


1.36/(10)® reyns 


As shown in Fig. 3.15, the gear axis whirl is essentially 
eliminated when the considered oil film bearing is used to 
support the shaft. The coupling effect at the mesh causes 
torsional oscillations without the existence of any torsional 
excitation. The magnitude of this motion is practically 


negligible in this case as can be seen in Fig. 3.16. 
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Fig. 3.13 Gear Whirl for a Stable condition with 
Unbalance on Rotor 1 (Case 10) 
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Fig. 3.14 Gear Whirl for a Stable Condition with 
Two In-phased Unbalances (Case 11) 
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Fig. 3.15 Gear Whirl for the system with Two 
Unbalanced Rotors Supported on Fluid 
Film Bearings (Case 12) 
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Fig. 3.16 Deflection of Gear Tooth for Case 12 
Ty = .0017, ry = .00125 
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un ce 

Lateral resonance of the two rotor system can be induced 
when the rotating frequency of the shaft is the same as the 
equivalent lateral natural frequency. In this case the 


following value is used. 


Resonant Respo) 
A resonant case is investigated for the following 


conditions. 


4, = +00003/M, inch 


Wy = 5000 rad/sec 


y = 5029 rad/sec 

The whirl motion of the gear axis in this case is plotted 
in Fig. 3.17 which illustrates the linear build up of the 
motion in this case with a frequency equal to the equivalent 
whirl frequency. 


ec! i 2 on ti! 


Another whirl motion is generated by using an additional 
unbalance on the other rotor in order to suppress the 
amplitude of oscillation due to the original unbalance. The 
additional excitation on the other rotor due to the same 


amount of unbalance with 0 degree phase shift can reduce the 
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deflection of gear tooth resulting from the first unbalance. 
The additional unbalance suppresses the high frequency 
oscillation due to the tooth contact induced by the first 


unbalance as shown in Fig. 3.18. 


fer 


the 

Fluid film bearings are used to support the shafts in 
this case. The parameters of the oil film bearing are selected 
in the same way as in the previous cases. The following 


parameter values are used as in the other unstable cases. 


c= .006" 


B= 1.36/(10)® reyns 


As shown in Fig. 3.19, the gear whirl is stabilized due 
to the positive damping provided by the selected oil film 
bearing. The deflection of gear tooth also shows acceptable 
magnitude (Fig. 3.20). 

As in section 3.1, the large amplitude of the oscillation 
due to an unbalance can be suppressed by the additional 
unbalance on the other rotor with a phase shift of 0 degree. 
Further reduction of the whirl motion can be achieved 
utilizing the damping provided by the proper selection of the 
oil film bearing. It is interesting to note that, as in case 
of the gear system with fluctuating torque inputs, the 
additional excitation has its more significant effect on 


suppressing the amplitude of high frequency oscillation caused 
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Fig. 3.17 Gear Whirl at Resonance beans 29) 
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Fig. 3.18 Gear Whirl due to Two Unbalances on Two 
Rotors at Resonant Frequency (Case 14) 


Ty = -00125, ry = .00125 


76 


-0001805 


0.000160- 


0.000140- 


0.0001204 


9.000100- 


0.0000804 


dg (inch) 


00000604 


0.000040- 


9000204 


900000--—,_, ; 
0000 0020 0040 
time (sec) 


T 1 
0060 0080 


Fig. 3.19 Gear Whirl subject to Two Unbalances with 
Appropriate Oil Film Bearings (Case 15) 
tT, = :00125, ry = .00125 
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Fig. 3.20 Deflection of Gear Tooth for case 15 
Ty = -00125, ry = .00125 
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by the tooth contact, and the damping due to oil film bearing 
works mainly on suppressing the build-ups. 

System with Subharmonic Resonance 
The instability of the two rotor system considered in 
this section is generated by inducing a subharmonic resonance 


situation as follows: 


This condition is selected based on the data given in 


Fig. 8(a) of the previous work by Seireg (22). 


6. sponse wit! Unb: 
The subharmonic resonant response is investigated for 


the following conditions. 


@, = .00003/M, inch 
@, = 10000 rad/sec 
¥y = $029 rad/sec 


As shown in Fig. 3.21, the subharmonic resonance also 


shows build-up of the oscillatory motion. 


ase The Effect of an Addit: Unbalance a 

As in the case of resonance, the large amplitude of the 
high frequency oscillation due to tooth contact excited by 
one unbalance can be suppressed by an additional unbalance 


excitation when the phase shift is 0 degree (Fig. 3.22). 
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Case 18: Effect of the Oi] Film Bearings 
The bearing parameters are selected as discussed 
previously. For this case, the following parameter values are 


selected. 


c= .006" 


# =1.36/(10)8 reyns 


As shown in Fig. 3.23, the gear whirl is stabilized, and 
the amplitude of deflection of gear tooth also attains an 
acceptable magnitude (Fig. 3.24). 

In section 3.1, it has been demonstrated that the 
amplitude of torsional vibration due to an existing 
fluctuating torque can be suppressed by an additional torque 
excitation and the appropriate selection of the oil film 
bearing parameters. In this section, it has been shown that 
the lateral response due to an unbalance can be neutralized 
also by an additional unbalance imposed on the system and the 


appropriate selection of the oil film bearing parameters. 
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Fig. 3.21 Gear Whirl at subharmonic Resonance (Case 16) 
Ty = 00063, ry = .00125 
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Pig. 3.22 Gear Whirl due to Two Unbalances on Two 
Rotors at Subharmonic Resonance (Case 17) 


1, = 00063, Ty = .00125 
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Fig. 3.23 Gear Whirl Subject to Two Unbalances 
under Subharmonic Resonance with 
Appropriate Oil Film Bearings (Case 18) 
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Fig. 3.24 Deflection of Gear Tooth for Case 18 
Ty = -00063, Ty = .00125 
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CHAPTER 4 
THE USE OF ONE TYPE OF INSTABILITY TO SUPPRESS ANOTHER TYPE 


The two mass pinion-gear system (Fig. 2.1) is considered 
here. The shafts are originally supported on ball bearings. 
In this chapter, an unbalance is used to suppress the 
oscillation due to a torsional excitation. 

It is conceptually possible to suppress a torsional 
oscillation by a lateral oscillation since both the lateral 
and the torsional vibrations of a geared system interact at 
the points of tooth contact during the mesh action. 

The examples discussed in this chapter are considered to 
illustrate the feasibility of utilizing the concept in the 
design of high speed gear systems where instabilities are 
expected to occur. 

It should be noted that, in the previous chapter, the 
additional input is of the same type as the original 
disturbance. This implies that the effect of the oil film 
bearings on the individual responses due to the two inputs is 
basically the same. Therefore, oil film bearings can be added 
if necessary after investigating the effect of the additional 
input. However, when a different type of input is imposed to 


the system, the effect of the oil film bearings on the two 


a4 


a5 
responses may not be the same. The system should therefore be 
analyzed prior to the introduction of the additional input to 


evaluate the effect of the oil film bearings on the response. 


Stable syste! 


In this case, if the amplitude of the response is 
allowable, then no additional treatment is necessary. In this 
section, however, the possibility of suppressing the response 
amplitude due to an existing torsional excitation on one of 
the gears by introducing an unbalance on the other gear is 
investigated. 
case s 

The stable response of the system (Fig. 4.1) due to an 
oscillatory torque input is evaluated for the same system as 
in the previous chapter (Fig. 2.1). 

The computation is undertaken for the following 


conditions. 


= Wy (1000) Ib-in 
0, = 3750 rad/sec 
= 3750 rad/sec 


Wyo = 5012 rad/sec 
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se Effect of the 

As discussed before, the bearing parameters (length, oil 
viscosity, and clearance) are selected to give the amount of 
damping necessary to reduce the amplitude of oscillation of 
the pinion-gear system due to the torque input. The following 


parameter values are selected. 


© = .005" 


u= 1.36/(10)8 reyns 


The tooth spring deflection in this case is shown in Fig. 
4.2. It can be seen that damping from oil film has eliminated 
the beat in the response as well as reduced the peak 
amplitude. 


Se: on of a "Tuned" whir’ ts 


on the System Response 

It is necessary in this case that the response due to 
the imposed unbalance have same frequency as the response due 
to the torsional excitation. 

In order for the response due to the unbalance to exhibit 
the same frequency as the torsional response, the rotating 
frequency of the shaft is selected to be the same as the 
torsional excitation frequency. Also, the amplitudes of the 
response due to the two different excitations should be of 


similar magnitude. Considering these factors, the following 
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parameters of the lateral vibrations are selected in an 


iterative way to obtain a tuned whirl. 


+00009/M, inch 


3750 rad/sec 


The tooth spring deflection for the system supported on 
the same oil film bearings due to the imposed unbalance 
excitation is shown in Fig. 4.3. 

This unbalance can be introduced in the system with 180 
degrees of phase shift which is found to provide the most 
reduction in the amplitude of oscillation due to the torsional 
input. As shown in Fig. 4.4, the amplitude of the tooth spring 
deflection is suppressed by the tuned whirl, and the gear 


whirl exhibits decay instead of build-up (Fig. 4.5). 


Resonance in this example is obtained by selecting the 
excitation frequency to be the same as the nominal natural 


frequency of the pinion-gear system. i.e., 
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Fig. 4.1 Deflection of Gear Tooth at the Mesh due 
to a Torsional Excitation at Rotor 1 
(Case 19) Te, = -0017, Th, = -00125 
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Fig. 4.2 Deflection of Gear Tooth at the Mesh due to 
a Torsional Excitation at Rotor 1 with 
the Effect of Oil Film Bearings (Case 20) 
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Fig. 4.3 Deflection of Gear Tooth at the Mesh due to 
an Unbalance Excitation at Rotor 3 with 
the Effect of Oil Film Bearings (Case 21) 


Ty = +0017, ry = .00125 


91 


0.00040- 
0.000304 
0.000204 
0.000104 


gs! Pt AF 
; 


-0.00010- 


dt (inch) 


-0.000204 


-0.000304 


N00 1 
0000 0020 0040 0060 0080 
time (sec) 


Fig. 4.4 Deflection of Gear Tooth at the Mesh due to 
the Two Excitations with Proper Phase Shift 
with the Effect of Oil Film Bearings (Case 21) 
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Fig. 4.5 Gear Whirl for Case 21 
= .00125, 
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cas ant Res: 
The resonant response of the system due to a torque 
excitation builds up linearly with time as would be expected. 
The response shown in Fig. 4.6 is obtained for the 


following conditions. 


T, = Wy (1000) 1b-in 
Voy = Ung = 5012 rad/sec 
0, = 5000 rad/sec 
It can be seen from the figure that the response build- 
up is at a frequency which is the same as the nominal natural 


frequency of the pinion-gear system. 


ase 23: The Effect of the Oil Film Bearings 


The effect of the damping in the oil film bearing is 
examined first. Bearing length, oil viscosity, and clearance 
of the bearing are in the same way as in the previous cases 
to provide positive damping. For this case, the following 


parameter values are selected. 


c= .006" 


w= 1.36/(10)® reyns 


As shown in Fig. 4.7, the build-up is suppressed by the 


positive damping provided by the oil film. 
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the ™ " whir! 


ns eS) 

The resonant response of the same system due to the 
unbalance can be induced to have the same form as that due to 
the fluctuating torque input when the torque exciting 
frequency is the same as the rotational frequency. The 


following inputs are used in the analysis. 


e3 = -000012/M, inch 


5000 rad/sec 


wy = 5029 rad/sec 


Wey = Ung = 5012 rad/sec 


The calculated tooth spring deflection for the system 
supported on the same oil film bearing with tuned whirl is 
shown in Fig. 4.8. 

The oscillation due to an unbalance can be introduced to 
the system with a 61.9 degrees of phase lag which is 
determined iteratively to suppress the resonance due to the 
torque input. As shown in Fig. 4.9, the amplitude of the tooth 
spring deflection is significantly reduced by the additional 
tuned input. Also, the gear whirl caused by oil film bearing 
and the unbalance excitation attains a stable behavior with 


a limit cycle (Fig. 4.10). 
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Fig. 4.6 Deflection of Gear Tooth at the Mesh due 
to a Torsional Excitation at Resonant 
Frequency at Rotor 1 (Case 22) 
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Fig. 4.7 Deflection of Gear Tooth at the Mesh due 
to a Torsional Excitation at Resonant 
Frequency at Rotor 1 with the Effect of 
Oil Film Bearings (Case 23) 
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Fig. 4.8 Deflection of Gear Tooth at the Mesh due 
to an Unbalance Excitation at Resonant 
Frequency at Rotor 3 with the Effect of 
Oil Film Bearings (Case 24) 
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Fig. 4.9 Deflection of Gear Tooth at the Mesh due to 
the Two Excitations at Resonant Frequency 
with Proper Phase Shift with the Effect of 
Oil Film Bearings (Case 24) 
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Fig. 4.10 Gear Whirl for Case 24 
Tex = 100125, Tho = -00125, 
ry%= .00125 | 7° 
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In this section, the same torsional instability which 
was discussed in section 3.1 is considered here to investigate 
the effect of introducing an unbalance on the rotor in order 
to suppress it. Again, the unstable situation, characterized 


by the exponential build-up, is expected when 


The same unstable response as in section 3.1. is obtained 
when the condition which satisfies Eq. 3.2 exists between 
excitation frequency of the torque, the nominal torsional 
natural frequency of the pinion-gear system, and the mesh 


frequency. The following input parameters are used. 


T, = W,(1000) 1b-in 
Wg = 2500 rad/sec. 
Ono = 5012 rad/sec 
0, = 10000 rad/sec 


As shown in Fig 4.11, the instability is characterized 


by the exponential build-up pattern. 
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Build-up frequency is same as natural frequency even 

though the excitation frequency is one half of that frequency, 

which is expected in the considered case of superharmonic 

resonance. 

Case 26: Effect of the Oil Film Bearing: 

As in previous cases, the effect of damping in the oil 
film bearings on the response is investigated first. The 
bearing parameters are selected based on the same scheme as 
discussed in sections 3.1 and 3.2. For this case, the 


following parameter values are used. 


c= .004" 


w= 1.36/(10)® reyns 


As shown in Fig. 4.12, the response of the system reaches 


a limit cycle due to the influence of oil film. 


election of the "Tui 
on the System Respons: 


‘The response of the system supported on the same oil film 
bearing due to the unbalance excitation (Fig. 4.13) is 
expected to have the same vibration frequency as the response 


due to torque input for the following conditions. 


3 = -00006/M, inch 


= 2500 rad/sec 


5029 rad/sec 
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The additional oscillation created by the unbalance 
excitation is introduced in the system with a 32.8 degrees of 
phase lag which is determined iteratively. The resultant 
tooth spring deflection is shown in Fig. 4.14, and Fig. 4.15 
shows the gear whirl caused by the unbalance. 

In sections 3.1, and 3.2, the additional inputs used to 
improve the system stability are of the same type as those 
causing the problems. However, in this chapter, it is shown 
that different type of excitation can be used to suppress the 
adverse effects of other disturbances. In other words, it is 
possible to introduce an appropriate amount of unbalance with 
the proper phase shift can be introduced in the system to 
suppress large oscillatory amplitudes arising from resonances 
or instabilities generated by torque excitations. 

The investigated cases show that the potentially 
destructive effect from a torque excitation can be avoided by 
“tuned" whirl and the proper design of oil film bearing 
without causing a large amount of lateral gear movement. All 
the reported results have been obtained without considering 
the effect of the inherent damping in the mechanical system 
itself. 

It can therefore be stated that it is possible in many 
situations to rely on the appropriate coupling of torsional 
and lateral instabilities to neutralize the adverse effects 
which can result from either of them if they would occur 


separately. 
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Fig. 4.11 Deflection of Gear Tooth at the Mesh due 
to a Torsional Excitation at Rotor 1 
for Unstable Condition (Case 25) 
Tey = +0025, Tho = -00125 
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Fig. 4.12 Deflection of Gear Tooth at the Mesh due 
to a Torsional Excitation at Superharmonic 
Resonant Frequency at Rotor 1 with the Effect 
of Oil Film Bearings (Case 26) 
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Fig. 4.13 Deflection of Gear Tooth at the Mesh due 


to an Unbalance Excitation at Superharmonic 
Resonant Frequency at Rotor 3 with the Effect 
of Oil Film Bearings (Case 27) 
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Fig. 4.14 Deflection of Gear Tooth at the Mesh due to 
the Two Excitations at Superharmonic Resonant 
Frequencies with Proper Phase shift with the 
Effect of Oil Film Bearings (Case 27) 
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Fig. 4.15 Gear Whirl for Case 27 
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CHAPTER 5 


APPLICATION TO GEARED SYSTEMS WITH TWO OVERHANG MASSES 


In this chapter, a whirl excitation due to an unbalanced 
overhang mass is introduced to a system to neutralize the 
effect of a torsional instability induced by a torque 
excitation on the other overhang mass. The objective is to 
investigate whether the same approach can be applied to 
complex system with more degrees of freedom than the two rotor 
system which has been considered so far. 

Basically, the system analyzed in this chapter (Fig. 2.2) 
is the same as the one treated in previous chapters except for 
the complexity (six more degrees of freedom) generated by 
introducing two additional rotors. The main distinction is 
that the effect of the excitation induced on the overhang 
rotors reaches the contact points on the teeth with some lag 
which depends on the dynamic parameters of the system. 

In order to demonstrate the feasibility of the approach, 
only the unstable case of exponential build-up of the response 
will be investigated. 

The following are the main system parameters used to 


investigate the system used in the reported study. 
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@ = 20 degree 

N= 4200 rpm 

0, = 2 N/ 60 rad/sec 

2 = 0, (GR) rad/sec 


0, = 2 (# of teeth of pinion) rad/sec 
D, = 2.5" 
Dy = 2.5" 


4, (bearing diameter) = D, + 2.¢ inch 
1, (bearing length) = 2.5" 
=u, <3" 
= Ly = by = bg = 16.2" 
+172 slug 
+172 slug 


x 
4 


4865 lb/in 
E = 30 (106) 1b/in? 
Poisson's ratio = .292 


Unstable Response due t: 
on _an Overhang Mass 


ase 


ation 


An unstable response (Fig. 5.1) with exponential build- 
up is obtained for the following conditions which causes a 
superharmonic resonance of the pinion-gear system supported 


originally by ball bearings. 


= 300 lb-in 
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Wyo = 316 rad/sec 


Way = 316 rad/sec 


In this case, the torque is applied on an overhang mass 
and transmitted to the pinion-gear system through the shaft. 


se 


1 Effect of the 0: 

Once again, the effect of the oil film bearings on the 
response is investigated before introducing a tuned whirl. 
The oil film bearing can be designed in the same way as in 
the previous cases. For this case the following parameter 


values are selected. 


© = .002" 


m= 1.36/(10)® reyns 


As shown in Fig. 5.2, the tooth spring deflection is 
stabilized to a limit cycle due to the effect of the positive 
damping provided by the oil film bearings. 


3 


Selectio) the 
‘on the System Response 


m and ot 


As in the previous cases, a whirl response with the same 
frequency (Fig. 5.3) is necessary to suppress the instability 
due to the torque. In this chapter the tuned response is 
obtained from the proper amount of eccentricity on the other 
overhang mass with a 150 degrees of phase lag which is 


obtained in an iterative way. 


qin 


e, = -08/M, inch 


, = 440 rad/sec 


The effect of the unbalance introduced in an overhang 
mass as well as the effect of the torque on the other overhang 
mass is transmitted with phase lags to the tooth mesh where 
the lateral and the torsional oscillations interact. 

As in the previous case of the two rotor geared system, 
the parameters are modified to match the frequencies of the 
system response to the torque input and the unbalance 
excitations. 

Using the two excitations along with the proper amount 
of the phase shift, the system shows very low response 
amplitude (Figs. 5.4 and 5.5) 

It can therefore be seen that the additional tuned whirl 
induced by the unbalanced overhang mass and the proper design 
of the oil film bearings can eliminate the potentially 
destructive instability due to the torque excitation applied 


on the other overhang mass. 
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Fig. 5.2 Deflection of Gear Tooth at the Mesh due 

to a Torsional Excitation at Rotor 2 

with the Effect of Oil Film Bearings 
02 


(Case 29) 1, 


ex 


i Pag 


102 


114 


0.02005 


0.01504 


0.01004 


-00504 


00005 


00504 


-0.01004 


~0,01504 


~0.0200 Lae eae tears roan aces OCT Cah ae 


000 .020 .040 ,060 .080 .100 .120 .140 .160 .180 
time (see) 


Fig. 5.3 Deflection of Gear Tooth at the Mesh due to 
an Unbalance Excitation at Rotor 2 with 
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Tex = +92, Tho = -02 


Fig. 


0200 
000 .020 


02005 


01504 


01004 


00504 


00004 


00504 


01004 


01504 


11s 


T 


T 
40 


060 .080 .100 
time (sec) 


T 
120 


5.4 Deflection of Gear Tooth at the Mesh due 
to the Two Excitations with Proper 
Phase Shift with the Effect of Oil Film 
Bearings (Case 30) 
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CHAPTER 6 
GUIDELINES FOR SYSTEM DESIGN 


Based on the results of the reported study, the following 
is a suggested procedure for guiding the designer in selecting 
the system parameters which insure that the response of the 


system remains within acceptable limits. 


Stable System with Torque tat: 


1. Calculate the system response for the given operating 
condition. If the response exhibits no build-up, and the 
vibration amplitude is acceptable, no more treatment will be 
necessary. 

2. If the system is stable and the amplitude is more than 
can be tolerated, the use of fluid film bearing with 
appropriately selected parameters can be investigated as a 
means of reducing the vibration to an acceptable level. This 
is accomplished by selecting the parameters of the fluid film 
bearings supporting the shafts in such a way as to introduce 
positive damping in the system. This is done by changing the 
pearing length, the radial clearance, and the viscosity 


iteratively until the oscillation is suppressed to the 
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allowable limit. This is illustrated by the results of the 
Cases 19, 20. The introduction of a fluid film bearing with 


ly = 2-5", o = 1.36/10® reyns suppressed 20 % of 


the amplitude of deflection of the gear tooth (from .0004" 
to .00032"), and eliminated the beat from the response as 
shown in Figs. 4.1 and 4.2. 

3. If necessary, further reduction of the system response 
can be achieved by imposing an additional "tuned" excitation 
on the other rotor of the system supported on the oil film 
bearing selected in the last step. If it is possible to 
introduce a torque excitation in the system, an additional 
torque with same magnitude and 180 degrees phase shift can 
practically eliminate the system oscillation. This is 
illustrated in Cases 2, 3. The introduction of the same amount 
of torque with 180 degrees of phase shift suppressed 99 % of 
the amplitude of deflection of gear tooth (from .0003" to 
+000003") as shown in Figs 4.2 and 3.3. In case additional 
torque excitation is not practical, the introduction of an 
unbalance on the other gear with a phase shift which is 
determined by matching the peak amplitude of the responses due 
to the individual excitations can further reduce the system 
response. This is illustrated in Cases 19, 21. The 
introduction of an unbalance excitation corresponding to an 
eccentricity with e, = .00053" with 180 degrees of phase shift 
suppressed 90 % of the amplitude of deflection of the gear 


tooth (from .0003" to .00003") as shown in Figs. 4.2, 4.4. 


ug 


1250} 


1. The reduction of the system response can be achieved 
by imposing an additional "tuned" excitation on the other 
rotor. If it is possible to introduce a torsional excitation 
in the system, an additional torque with same magnitude and 
180 degrees phase shift can reduce the system response 
drastically. This is illustrated in cases 4, 5. The 
introduction of the same amount of torque with 180 degrees of 
phase shift eliminated 99 % of the amplitude of deflection of 
the gear tooth (from .0015" to .000002") at t =.007 sec. as 
shown in Figs. 3.5 and 3.6. 

2. The introduction of a fluid film bearing can completely 
eliminate any residual effect. The procedure of the selection 
of the parameters of the oil film bearing is the same as 
discussed in the stable case. This is illustrated in Cases 5, 


6. As shown in Figs. 3.6 and 3.7, the oil film bearing with 


ly = 2.5" +006", » = 1.36/108 reyns reduced the amplitude 
of deflection of the gear tooth to 0.02 % of the original 
amplitude (from .000002" to .0000003"), and eliminated the 
build-up. 

In case it is not possible to impose an additional torque 
on the system, the use of an unbalance in the other gear with 
a phase shift which is determined by matching the peak 


amplitude of the responses due to individual excitations can 


reduce the system response. In this case, the effect of oil 
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film bearings on the system response subject to each 
excitation is different. Accordingly, the individual response 
of the system supported on the selected bearings due to each 
excitation is obtained first to decide the appropriate amount 
of phase shift and unbalance. 

The magnitude of unbalance is determined by matching the 
peak amplitude and the amount of phase shift is selected in 
such a way as to make the peak responses out of phase with 
each other. This is illustrated in Cases 22, 23, and 24. The 
response of the system supported on the appropriately selected 
bearings due to torsional and unbalance excitations 
respectively is shown in Figs. 4.7 and 4.8. The introduction 
of an unbalance excitation corresponding to an eccentricity 


@3 = -000012" with 61.9 degrees of phase lag along with 


properly designed bearings (1, = 2.5", c= .006", 4 = .5/10%) 
suppressed 95 % of the amplitude of deflection of gear tooth 
(from .002" to .0001"), and eliminated the build-up as shown 


in Figs. 4.6 and 4.9. 


stable System with Torque Excitation 


1. The stabilization of the system response with 
exponential build-up can also be achieved by imposing an 
additional "tuned" excitation on the other rotor. If it is 
possible to introduce a torsional excitation in the system, 


an additional torque with the same magnitude and 180 degrees 
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phase shift can reduce the system response drastically. This 
is illustrated in Cases 7, 8. The introduction of the same 
amount of torque with 180 degrees of phase shift eliminated 


99.9 % of the amplitude of deflection of the gear tooth (from 


+005" to .000003") at t =.007 sec. as shown in Figs. 3.9 and 
3.210. 

2. The introduction of a fluid film bearing can further 
reduce the response. Again the parameters of the oil film 
bearings are selected in such a way to provide the system with 


positive damping as discussed in the previous cases. This is 


illustrated in cases 8, 


As shown in Figs. 3.10 and 3.11, 


the oil film bearing with 1, = 2.5", c = .004", w = 1.36/10% 
reyns can further suppress the deflection of the gear tooth. 

once again, in case it is not possible to impose an 
additional torque on the system, the introduction of an 
unbalance on the other rotor can reduce the system response. 
The procedure for selection of the magnitude of unbalance and 
the amount of phase shift is the same as in the resonant case. 
This is illustrated in cases 25, 26, and 27. The response of 
the system supported on the appropriately selected bearings 
due to torsional and unbalance excitations respectively is 
shown in Figs. 4.7 and 4.8. The introduction of an unbalance 
excitation equivalent to an eccentricity e, = .00006" with a 
32.8 degrees of phase lag along with properly designed 


= 004", » = 1.36/10® reyns) suppressed 


bearings (1, = 2.5" 


95 % of the amplitude of deflection of gear tooth (from .004" 
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to .0002"), and eliminated the build-up as shown in Figs. 4.11 


and 4.14. 


Stable System with Unbalance Excitation 


1. Calculate the system response for the given operating 
condition. If the response exhibits no build-up, and the 
vibration amplitude is acceptable, no more treatment will be 
necessary. 

2. If the system is stable and the amplitude is more than 
can be tolerated, the reduction of the system response can be 
achieved by imposing an additional "tuned" excitation on the 
other rotor. An additional unbalance excitation with the same 
magnitude and 0 degree phase shift can further reduce the 
system response. This is illustrated in Cases 10, 11. The 
introduction of the same amount of unbalance with 0 degree of 
phase shift suppressed approximately 90 % of the amplitude of 
the high frequency gear whirl (from .0003" to .00003") as 
shown in Figs 3.13 and 3.14. 

3. The introduction of fluid film bearings selected in an 
iterative way to supply the proper amount of damping to the 
system can further stabilize the system response. This is 
illustrated in Cases 11, 12. As shown in Figs. 3.14 and 3.15, 
the oil film bearing with 1, = 2.5", c = .005", u = 1.36/10° 
reyns eliminated 95 % of the fluctuation in gear whirl (from 


+0004" to .000002"). 
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soni 0 Unbalance 


1. The reduction of the system response can be achieved 
by imposing an additional "tuned" unbalance excitation on the 
other rotor. An additional unbalance excitation with same 
magnitude and 0 degree phase shift can reduce the system 
response drastically. This is illustrated in Cases 13, 14. The 
introduction of the same amount of unbalance with 0 degree of 
phase shift eliminated approximately 90 % of the amplitude of 
the high frequency oscillation of gear whirl (from .0012" to 
+0001") as shown in Figs. 3.17 and 3.18. 

2. The introduction of fluid film bearings selected as 
discussed previously can stabilize the system response. This 


is illustrated in Cases 14, 15. As shown in Figs. 3.18 and 


3.19, the oil film bearings with 1, = 2.5", c = .006", w= 
1.36/10® reyns eliminated the build-up, and stabilized the 


response at .00035". 


Subharmonic Resonance due to Unbalance Excitation 


1, The reduction of the system response under subharmonic 


resonance can also be achieved by imposing an additional 


tuned" unbalance excitation on the other rotor. An additional 
unbalance excitation with the same magnitude and 0 degree 
phase shift can reduce the system response drastically. This 


is illustrated in cases 16, 17. The introduction of the same 
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amount of unbalance with 0 degree of phase shift reduced the 
amplitude of the high frequency oscillation of gear whirl from 
+005" to .0006" as shown in Figs. 3.21 and 3.22. 


2. The introduction of fluid film bearings which is 


selected in an iterative way to supply positive damping to 
the system can enhance the stability and further reduce the 
response. This is illustrated in Cases 17, 18. As shown in 


Figs. 3.22 and 3.23, the oil film bearings with 1, = 2. 


.004", 4 = 1.36/10® reyns caused the amplitude of gear whirl 


to decay. 


CHAPTER 7 
CONCLUSIONS AND RECOMMENDATIONS 


This study demonstrates the feasibility of utilizing 
lateral and torsional instability properties in multi-rotor 
geared system to neutralize the potentially destructive 
response if they occur separately. The use of oil film 
bearings to stabilize and suppress the amplitude of response 
has been illustrated by case studies. A procedure is given to 
guide the designer in eliminating the instability potential 
in the design stage by utilizing other instabilities and 
damping in the bearings. 

The findings of the reported investigation can be 
illustrated by the case studies summarized in Table 6.1. 

Many cases are considered to illustrate the different 
possibilities of suppressing the amplitudes of system 
response. In Chapter 3, additional excitations of the same 
type as the existing disturbances are imposed on the system 
to suppress the response. 

As shown in section 3.1, the response of the system due 
to a fluctuating torque excitation on a rotor can be 
suppressed by introducing additional fluctuating torque input 


with 180 degrees phase shift in the system at the other rotor. 
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The response of the system can be further suppressed by the 
damping provided by appropriate selection of the parameters 
of oil film bearings. The combined result of imposing an 
additional input and the appropriate selection of the oil film 
bearing parameters can cause the system response to exhibit 


a decay or a limit cycle under conditions of stable, resonant, 


and unstable system behavior as shown in Cases 3, 6, and 9 
respectively. 

The possibility of suppressing the system response due 
to an unbalance on one rotor by an additional unbalance on 
the other rotor is investigated in section 3.2. An additional 
unbalance excitation with 0 degree phase shift along with the 
properly selected oil film bearings can reduce the systen 
response and eliminate the possible instabilities which may 
occur due to an existing unbalance as shown in Cases 12, 15, 
and 18 respectively. 

The discussions in Chapters 4 and 5 show that it is 
possible to suppress the effect of an existing disturbance by 
imposing a “tuned” additional excitation of a different type. 
It is shown in Chapter 4 that the system response due to an 
existing torsional excitation on a rotor can be suppressed by 
an appropriate amount of unbalance imposed on the other rotor 
along with the appropriate selection of the parameters of the 
oil film bearings. The amount of unbalance and the phase shift 
for the unbalance excitation are determined iteratively. An 


additional tuned excitation of a different type along with the 
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appropriate design of the oil film bearings can stabilize the 
system and reduce the response amplitude as shown in Cases 21, 
24, and 27. 

In Chapter 5, it is shown that the same approach can be 
applied to complex systems with more degrees of freedom. To 
demonstrate the feasibility of the approach, a worst case 
condition is investigated (Case 28). The instability due to 
an existing torsional excitation on an overhang mass is 
neutralized by imposing an unbalance on the other overhang 
mass along with appropriate selection of the parameters of 
oil film bearing as shown in Case 30. 

It can be concluded from this study that "tuned" 
instabilities can be induced in geared systems to provide safe 
operation at critical speeds where resonances and 
instabilities may occur. It is therefore possible to design 
a multi degree of freedom geared system where the build-up of 
the forces at the mesh can be avoided by appropriate selection 
of the system parameters. This is a practical approach since 
both the lateral and the torsional oscillations in gear 
systems interact at the points of tooth contact. 

The suggested procedure in Chapter 6 can be used to guide 
the designer in selecting the system parameters which insure 
the stability of the system under the expected critical 


conditions. 
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Table 7.1 Effect of the "Tuned" Input and Oil Film Bearings 


Excitation 


Initial 
response 


Induced 
change 


Final 
response 


Limit cycle 


limit cycle 


| a=.0004" | a=az-o7" | 
(Fig.3.1) (Figs.3.3,3.4) 
Linear decay 

| build-up c=.006 

a=.25 (t) BL. 36/10" 

(Fig. 3.5) (Figs. 3.7, 3.8) 
1,=1500 Expon. T=-7, limit cycle 
,,=2500 | build-up 004 | a,=38-06" I 
(w,,=5012) d=e*-1 +36/10° 

(Fig. 3.9) (Figs. 3.11, 12) 
Me,=3E-05 Oscillating | e,=e, constant 
o=3750 | d.=.0004" c=.005 | a,=.0004' 
(,=5029) (Fig.3.13) u=1.36/10° | (Figs. 3.15, 16) 
Me=3E-05 | Linear constant 

| build-up | dj=-o004" = | 
slope =.40 

(Fig.3.17) (Figs. 3.19, 20) 

Linear constant 

| build-up | a,=.0003" i 
slope=.02 

(Fig. 3.21) (Figs. 3.23, 24) 


Limit cycle 


limit cycle 


| d=.0004" | a=4E-04" 
(Fig. 4.2) (Figs. 4.4, 5) 
Linear limit cycle 

| build-up | a=1E-04" 

| a=.25 (t) H=1.36/10° | 
(Fig. 4.6) phase -.34n | (Figs. 4.9, 10) 
Expon. limit cycle 

| build-up | E04 

| d=e*t-1 | 


| (Pig.4.12) 


(Figs. 4.14, 


if 
15) 
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Future Work 

The following are recommendations for future work before 
the proposed approach can be used with certainty in practical 
applications. 

Accurate evaluation of the film parameters is central to 
the reliable modelling of the whirl motion and is consequently 
an essential element of the developed procedure. The behavior 
of the fluid film is assumed to be isoviscous in the reported 
study. Seireg et al (3, 19, 45] show that the stiffness and 
damping characteristics can be significantly different when 
the thermohydrodynamic behavior of the bearing is considered. 
They suggested an empirical numerical procedure for computing 
the transient thermohydrodynamic restoring and damping forces 
generated in the film. This can have considerable effect on 
the stability threshold as well as the amplitude and frequency 
of the whirl motion. 

Experimental verification of the suggested design 
procedure is necessary to insure the validity of the approach. 
This is particularly important in view of the dependency of 
the analysis on the prediction of the behavior of thin oil 
films (which are rheological in nature), the out-of-roundness 
and the surface roughness of the journal and sleeve as well 
as the thermal and chemical characteristics of the lubricating 
oil. The thin chemical and absorbed layers generated on the 
bearing surface in the course of its useful life can also have 


a significant influence on the behavior of the bearing. 


APPENDIX 


List of Programs 


eacananaa 


0 


e990 
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TWELVE DEGREES OF FREEDOM GEARED SYSTEM 


* 
* 
++ MAIN PROGRAM. --- * 
FOR THE CORNELL SUPER COMPUTER * 

* 


DOUBLE PRECISION 


IMPLICIT REAL*8 (A-H) 
IMPLICIT REAL*8 (0-2) 


REAL*8 ¥1(60010) , X2(60010) ,x3(60010) ,x4(60010) 
REAL*8 ¥1(60010) , ¥2(60010) ,¥3(60010) , ¥4(60010) 
REAL*S ‘TH1(60010) ,TH2(60010) , TH3(60010) ,TH4(60010) 
REALS ‘DP(61000) ,DG(61000) , DTS(61000) 


REAL 11,12,J1,J2,J3,J4,NU,121,122,KT1,KT2, INCRE 
REAL L1,12,13,14,15,16,K0,K, JEFF 


COMMON WE1,WE2,WE3,WE4,D1,D2,1,L2,14,15,GR 
COMMON BL1,BL2,BD1,BD2,C,W1,W2,N 


OPEN(UNIT=30, STATUS="NEW" , FILE="MAIN. DAT" ) 


GIVEN PARAMETERS- 
D'S FOR SHAFTS, R'S FOR ROTORS 


2T=.0 
CONST=.0 
CONV=180/3..14159 
CONV1=386. 
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PSI1=0. /CONV 
PSI2=0. /CONV 
PSI3=0. /CONV 
PSI4=170. /CONV 
‘X12=0. /CONV 

X14=0. /CONV 


WI=NK3.14159#2/60 
W2=GRAW1 

WEI=3. 14159*RIM*2HWIX, 282 
WE2=WEL 
WES=3.14159*RI*#2*WI*. 262 
WEO=WE3 

‘XM1=WE1/CONVI 
XM2=WE2/CONVL 
XM3=WE3/CONVI 
XMG=WE4/CONVI 

E1=.0/xM1 
E3=.0/XM3 
E2=.0/XM2 


I1=3, 14159*D1**4/64. 
12=3.14159*D2**6/64. 
TA=KMIAR D2 /2. 
J2=XM2#R22/2. 
JS=XMSHRIN2 / AGRA, 
Tu=KMGeRUsr2 /24GRIHD, 
NU=.292 
G=E/(2#(14NU)) 


(C--- INFLUENCE COEFFICTENT-~ 


2463 /(6*EAT1) 

DUAL2*#2) /(4*E*11) 

eH 2#(L14 2412) /(3*E*11) 
3 / (6*E*I2) 
LA*LS##2)/(4*E*12) 

Late De (Lt 2#LS) /(3*E*12) 


CONSTANT RELATED TO THE POSITION OF BEARING~ 


133 


BLI=(L1*#24(L14L2#2)##2) /(L14L2)**2 
2RL HD) /(LIFL2)*#2 


12 
B33=(L4##2+(IA4LS*2)*42) /(LA4L5 )**2 


B3G=(2ALSHH2) / (14415 482 
BG4=B34 


© CONSTANT RELATED TO THE MAG. OF EIGENVALUES 


PU=PAIXG 
PS=PATY] 
P6=PAIY2 
P7=PAIYS 


JJIN DECIDES SAMPLING INTERVAL OF OUTPUT FILE. 


c 
KK DECIDES MULT. OF CALC. TIME BY ADDING DO L00P(7000). 


¢ 


‘J4J=20000 


© CCR IS ACCURACY CRITERION FOR ITERATION 
‘ACCR=.0001 


DTMIN=. 05%3, 14159/180/(P#IMAX) 
INCRE=PATMAX#DTMIN#180. /3. 14159 
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DTJJ=KK*IMAX*DTMIN 


DT1=DTMIN*IMAX/TX1 
DT2=DTMIN*IMAX/1X2 
DT3=DTMIN*IMAX/ IX3 
DTG=DTMIN*IMAX/ IX4 
DIS=DTHIN*IMAX/TY1 
DT6=DTMIN*IMAX/ 1Y2 
DT7=DTMIN*IMAX/ IY3 
DT=DTMIN*IMAX/1¥4 
DT9=DTMIN®IMAX/ ITH 
DT10=DTMIN*IMAX/ITH2 
-DTMIN*IMAX/ITH3, 


© ADT GIVES THE INCREMENTAL ANGLE ON EACH PHASE PLANE 


K0=4865 
-KO*RIM#2#COS( PHI )*#2 
5 
JEFF=(J1*J3)/(J14I3) 
WNO=(K/JEFF)*. 5 


TZ1=3. 14159*D 1444/32 
122=3.16159*D2**4/32 
KTI=IZ1*G/(L14L2) 
KT2=1224G/(L44L5)*GR**2 


© INIT. CONDITIONS 
x1(1)=.0 


v2(1)=.0 
yoT=.0 


¥3(1)=.0 
y3T=.0 
Dy3=.0 
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LCL) **2+YI (1 )eAD YAH, S 
:3(1)#H24YS(1) #2 JH, 5 
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CALL SUBROUTINE FOR K'S & C'S DUE TO BEARING- 


CALL BRNG(S1,$2,XXK1,XXK2, XYK1, XYK2, YXK1, YXK2, YYK, YYK2,, 
* ‘CXX1, CXX2, CXY1, CXY2, CYK1, CYK2,CY¥1,CYY2) 


‘DO LOOP FOR INTEGRATION STEPS BY PHASE PLANE METHOD- 


THIS DO LOOP(7000) IS MADE JUST TO SEE THE DISPLACEMENTS OF THE 
SYSTEM LONGER W/O VIOLATING THE VIRTUAL MEMORY SIZE. 

X(JJ+1) IS THE DISPLACEMENT AFTER KK TIMES OF CALCULATION 
(EACH WITH INCREMENTAL ANGLE PADT) SINCE X(JJ). 


eane 


DO 7000 JJ: 


IIT 


DO 6000 II=1 


KK 
DO 6100 NM=1, 1X1 


SM1=KO#(1-A*COS(WH#T1)) 
DMI=2T#2( SMI*XMS#XMI/ (XMI4XML) J, 5 


CALL SIGN(DX1,SGN) 


DIX1=(XIP*(COS( PHI) **2*SM14XXKI*B1141 /A11)+ 
DX1(COS( PHI )**2*DM14CXX1B11)+ 
‘YIT#(COS(PHI)*SIN(PHI) *SM5+XYK1*B11)+ 
DY1*(COS(PHI)*SIN( PHI) *DM5-+CXY1*B11)+ 
‘THITRI*COS( PHI) #*2*SN9+ 
DTHI*R1*COS(PHI)**2*DH9- 
X2T*(XXKIMB12+1/A12)~ 
DX2*CXK1#B12- 

Y2T*XYKI*BI2~ 
DY2*OXY1*B12- 
XST*COS(PHI)*A2eSHS~ 
DX3*COS(PHI)**28DH3~ 
YST®COS(PHI)*SIN(PHI)*S47- 
DY3*COS(PHI)*SIN( PHI) *DM7-+ 
‘TH3T*R34COS( PHI) **2*SH1 1+ 
DTH3*RI*COS(PHI)**24DM11 
++SGNAXMIAW1*2*ABS(XIT)*CONST 
~XMIMELAY1442#C0S (WI*T14PS11) 
~A22/A12*XM2*ED#W1**2*COS(WIATI4PSIZ)) 
/CMIEPL##2)-X1T 


eee eee eee eee 
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‘CXAT+DTX1) 2+ (DX1/P1)*#2)**, 5 

PARXL*SIN(P1*DTI/2) 

CALL PSIS(X1T,DX1,DTX1,PSX1,P1) 

CALL PHASE(X1T, DXi, DTX1, PSXi, ANG, SGN1, SGN2, P1,DT1) 
T 


X1T=XIT#SGN1*XL1*SIN(ANG) 
DX1=(DX1/P1+SGN2*XL1*COS( ANG) )*P1 


c ITERATION 


6101 XB1=(X1T#XITC)/2 
DxB1=(DX14DX1C) /2 


CALL SIGN(DXB1, SGN) 


DITXBI=(XB1#(COS (PHI )*#2SM14XXK1*B1141/A11)+ 
DXBI*(COS( PHI) #+2*DMIHCXXI*B11)+ 
YBI1*(COS( PHI) #*SIN(PHI)*SH5+XYK1#B11)+ 
DYB1*(COS( PHT) *SIN(PHI)*DMS+CXY1#B11)+ 
‘THB1*R1*COS (PHI) *#26S49+ 
DIHBI*RIACOS(PHI)*#2*DM9- 
XB2*(XXK1*B124+1/A12)~ 
DXB2*CXX1*B12- 

‘YB2*XYKI*B12- 
DYB2*CXY1*B12- 
XB3*COS(PHI)**24SM3- 
DXB3*COS(PHT)**24DH3~ 
‘YB3*COS(PHI)#*SIN(PHI)*S47~ 
DYB3*COS( PHT) #*SIN(PHI)*DM7+ 
‘THB3*R3*COS (PHI) *#2*SHLL+ 
DTHB3*R3*COS(PHI)**24DH11 
+SGNXM LW *ABS (XB1)*CONST 
~XMLELAW1*2*COS(WI#TI+PSI1) 
7A22/A124KM2*E2*W1%#2*COS (WIATI4PSI2)) 
/(SMIAPI¥A2) “XBL 


eee eee ee eee ees 


RXB1=((X1TCHDTXBL) H+ (DXIC/P1) #42) #5 

XLB1=24RXBI*SIN(PIADT1/2) 

CALL PSIS(X1TC, DX1C, DTXB1, PSXB1, P1) 

CALL PHASE(X1TC, DX1C, DTXBi, PSXB1, ANG, SGN1, SGN2,P1,DT1) 
XBI=XITC#SGNI#XLB14S1N( ANG) 
DXB1=(DXIG/P1+SGN2*KLB1*COS(ANG) )*P1 

DIX1=(ABS(X1T) ~ABS(XB1) )**24(ABS(DX1) 
[PL -ABS(DXB1)/P1) #2 
TFCDIXL. GT. ACCR#(X1T#2+(DX1/P1)*#2) )THEN 

XIT=XBI 

DX1=DXB1 

GO TO 6101 


6102 
6100 


eee eres: 


ELSE 
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GO TO 6102 


ENDIF 


Ti=t40T1 


cr 


Do 6200 


CALL SI 


;ONTINUE, 
mete, 1X2 


GN(DX2, SGN) 


DTX2= (-XYT#(XXK1*B1241/A12)- 


ITERATION 


6201 


eee e res 


DXI*CKXI*B12~ 

YITAXYKI*B12~ 

DYI*CXYI*B124+ 

X27T*(XXKI*B22+1/A22)+ 

DX2eCxXI¥B22+ 

YOTRXYKI*B22+ 

Dy2*OXY1#B22 
-FSGNAAHZAW 12 ABS (X2T)*CONST 
~KH2HE2HW14424COS(WI*TZ4PSI2) 
~ALL/A124XM IE 14W14424COS(WIAT2+PST1)) 

/(8MQ*PIHKD) -K2T 


(X2THDTX2)*424(DN2/P2)¢AZ) IHS 
‘XL2=24RX2#SIN(P2*DT2/2) 
CALL PSIS(X2T, DX2,DTX2, PSX2,P2) 
CALL PHASE(X2T, DX2,DTX2, PSX2, ANG, SGN1, SGN2,P2,DT2) 
X2TO=X2T 
DX2C=Dx2 
X2T=X2T+SGN1*XL2*SIN(ANG) 
DX2=(DX2/P2+SGN2#XL2*COS(ANG))*P2 


XB2=(X2T#X2TC) /2 


DXB: 


CALL 


DIxp2: 


‘DX24DX2) /2 
SIGN(DXB2, SGN) 


(XBL (XXK14B1241/A12)~ 

DXBIRCXK1*B12~ 

YBIAXYKI*B12~ 

DYBIACXY1*B12+ 
XB2*(XXKI*B2241/A22)+ 
DXB2ACXXI*B22+ 

YB2AXYKINB22+ 

pyB2Acxy1*B22 
+¥SGN®XM2#W1#2:ABS(XB2)*CONST 
= XM2#EZ#W1*2*COS (WI*T24PSI2) 
~ALL/A124XH IKE AW 14424COS (WIAT2+PSI1)) 
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* / OMe 2mD) -XB2 


(X2TCHDTXBZ)*#2+(DK2C/P2)*42)**. 5 

XLB2=2#RXB2*SIN(P2*DT2/2) 

CALL PSIS(X2TC, DX2C, DTXB2, PSXB2, P2) 

CALL PHASE(X2TC, DX2C, DIXB2, PSKB2, ANG, SGN1, SGN2,P2,DT2) 
‘XB2=X2TC+SGN1*XLB2*SIN(ANG) 
DXB2=(DX2C/P2+SGN2*KLB2*COS (ANG) )#P2 

DLX2=(ABS(X2T) -ABS(XB2) )*#2+(ABS(DX2) 

* — /P2-ABS(DXB2)/P2)**2 
TF(DLX2. GT. ACGR#(X2T#*2+(DX2/P2)#2)) THEN 
XQT=XB2 
DX2=DXB2 
Go TO 6201 
ELSE 
Go TO 6202 
ENDIF 


6202 ‘T2=7240T2 
6200 CONTINUE 


DO 6300 MMMM=1, 1X3 


'SM3=KO#(1-A¥COS(WH¥TS) ) 
DM3=2T*22( SH3*XMKMI / (XM34XH1) 4,5 


CALL SIGN(DX3,SGN) 


DIX: 


“XI TACOS(PHI)*#2SH1- 
DxT#COS(PHT)##2#DM1~ 
YXD*COS(PHI)#*SIN(PHI)*SHS~ 
DY1*COS(PHI)*SIN(PHI)ADMS- 
‘THYTARIACOS(PHI)**2*SH9~ 
DIWI*RI+COS(PHT)**2*DHIE 
XST#(COS(PHL)*#2*SH3+XXK2*B3341 /A33)+ 
DX3#(COS (PHT) #**2*DM3+CXX2*B33)+ 
‘Y3T#(COS(PHI)*SIN(PHI)*SH7+XYK2#B93)+ 
DY3*(COS(PHT)*SIN(PHI)*DH74CXY2#B33) ~ 
THST*RI*COS (PHI) **24SHI1~ 
DTH3*RI4COS (PHI) **2ADHI1- 
XGT¥(XXK2*B34#1/A34)~ 

DX4ACKX2*B34~ 

YOT#RYR2#B34~ 

DY6ACKY2#B34 
-+SGNHXUSHW24ORABS(X3T)*CONST 
-XUS#ESAW2##24COS(W2AT34PSI3) 

Ab / ASGANMGAEGARD#HZACOS(W2*TS#PSI4)) 
/M3*P3402)-X3T 


ee eee eee ee eee es 


RX3=( (X3THDTXS)*#2+(DX3/P3)#H2)#, 5 
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XL3=2#RX3#SIN(P3*DT3/2) 
CALL PSIS(X3T,DX3, DTX3, PSx3,P3) 
CALL PHASE(X3T, DX3, DIX3, PSX3, ANG, SGN1, SGN2,P3,DT3) 
X3TO=X3T 
DX3C=Dx3 
XST=X3T+SGNI*XL3*SIN(ANG) 
DX3=(DX3/P3+SGN2*XL3*COS (ANG) )P3 


ITERATION 


6301 


DXB3=(DX3+DX3C)/2 
CALL SIGN(DXB3, SGN) 


DITXB3=(-XB1*COS( PHT) **2*5M1- 
DXB1*COS(PHI)*#2*DM1- 
YBI*COS(PHI)#SIN( PHI) #SHS~ 
‘DYB1*COS( PHI) *SIN(PHI)*DM5~ 
‘THBI*R1*COS(PHI)**2*SH9~ 
DTHB1*RI*COS(PHT)*#2*DHI+ 
‘XB3#(COS( PHI) *#2#SM3+XXK2*B3341 /A33)+ 
DXB3*(COS( PHI )**2*DM3+CXX2#B33)+ 
‘YB3*(COS( PHI) #SIN( PHI )*SH7-+XYK2*B33)+ 
DyB3¥(COS(PHI)*SIN(PHI)*DM7+CXY2*833) ~ 
‘THBS*R3*COS( PHT) **2SM11- 
DTHB3*R3*COS( PHI) **2*DM11- 

BQH (XXK2*B3441/A34) ~ 
DXBG#CKX2*B34~ 
YBOAXYR2*B34~ 
DYBGHCKY2*B34 
++SGN# 013 #W/24#2ABS( XB3)*CONST 
-KMS*E3*W2#*2C0S (W2*T3+PSI3) 
AGG /A3G#XMGHEGHW2*#2*COS(W2ATS+PSI4)) 
/(MBAP3##2) -XB3 


eee eee ee eee eee ee 


S 


(X3TCHDTXAI)#A2+(DKIC/PI) HD), 5 

XLB3=2*RXB3*SIN( P34#DT3/2) 

CALL PSIS(X3TC, DX3C, DTXB3, PSXB3,P3) 

CALL PHASE(X3TC, DX3C, DTXB3, PSXB3, ANG, SGN1, SGN2, P3,DT3) 
XBI=X3TC+SGNI*XLBS#SINCANG) 
DXB3=(DX3C/P3+SGN2*XLB3*COS(ANG))*P3 

DLX3=(ABS(X3T) -ABS(XB3))*#24(ABS(DX3) 

* /P3-ABS(DXB3)/P3)##2 

TF(DLX3, GT. ACCR*(X3T##24+(DX3/P3)*#2) )THEN 

X9T=XB3 


Go TO 6301 
ELSE 

GO TO 6302 
ENDIF 


aan 


6302 T3=T34DTS 
6300 CONTINUE 


DO 6400 MMMMEH=1, 1x6 
CALL SIGN(DX4, SGN) 


DIX4=(-X3T*(XXK2*B3441/A34)- 

DX3*CKX2*B34- 

YSTAKYK2#B34~ 

DYSACKY2*B34+ 

KOT (XXK2#B 4441 /AG4)+ 

DXGACXKD*RGG+ 

YOTAXYK2*BG4+ 

DY4ACKY2*BGG 
-+SGNAXMGAW2*424ABS(XOT)*CONST 
-KMGHEGHH2#%2*COS(W2ATG+PST4) 
-A33/A34#XN3*E3*W2**2*COS (W2*T44PSI3)) 

/CRMGHPGHE2) -KOT 


Reese eee 


RXG=((RUTHDTXG ) #24 (DK4/P4) #2), 5 
XLA=24RXGASIN(PGADTG/2) 
CALL PSIS(X4T, DX4,, DIX4, PSX4,P4) 
CALL PHASE(X4T, DX4, DIX4, PSX4, ANG, SGN1, SGN2, P4,DT4) 
XOTO=X6T 
DX4C=DX4 
XGT=X6T#SGNI*XL4*SIN(ANG) 
‘DX4/P4+SGN2*XL4*COS( ANG) )*P4 


ITERATION 


6401 
DXB4=(DXG#DX4C) /2 


CALL SIGN(DXB4, SGN) 


DTXB4=( -KB3*(XXK2*B3441/A34)- 
DXBSACKX2*B34- 
YBS*XYK2*B34~ 
DYBSACKY2*B34+ 
XBO*(XKK2*B4G41/A66)+ 
DxBOACKKZABAGE 
YBO#XYR2*B44+ 
DYBUACKY2#BAG 
4+SGNFAHGANDHE2RABS (BA) CONST 
“XMGAEGAW2##2*COS (W2*TH#PSTS) 
~A33/A34*XUSFES*W2**2*COS (W2*TU+PSI3)) 
MPI) -XBG 


eee ees 


RXBG=( (KGTCHDTABS )##2+(DK4C/PA)*RZ)*E. 5 
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XLB4=2*RXB4KSIN(P4*DTS/2) 

CALL PSIS(X4TC, DX4C, DTXB4 , PSXB4, P4) 

CALL PHASE(X4TC, DX4C, DTXB4, PSKB4, ANG, SGN1,SGN2, P4,DT4) 
‘XB4=KGTC+SGNI*XLBG*SINCANG) 
DXB4=(DX4C/P4#SGN2*XLB4*COS( ANG) )*P4 

DLX4=(ABS(X4T) ~ABS(XB4) )*#2+(ABS(DX4) 

* _ /P4-ABS(DXB4)/P4)*42 
IR(DLX4. GT. ACCR#(X4T#'2+(DX6/P4) #2) )THEN 
XGT=XB4 
DX4=DXB4 
GO TO 6401 
ELSE 
GO TO 6402 
ENDIF 


6402 TH=T4HDTS 
6400 CONTINUE 


DO 6500 LL=1,1¥1 


SMS=KO#( 1-A*COS(WH*TS )) 
DMS=27*2%( SHS*XM3*XM1 / (XM34XML) 4. 5 


CALL SIGN(DY1,SGN) 


DIY1=(X1T*(COS(PHI)#SIN(PHI)*SM14YXKI8B11)+ 
DX1*(COS( PHI) #*SIN(PHI)*DM14CYX19B11)+ 
YET#(SIN(PHI)*#2#SM5-+YYKI*B1141 /A11)+ 
DY1*(SIN(PHI)#2#DM5+CYY18B11)+ 
THITARIMSIN( PHT) #42#SM94+ 
DTWI*RI*SIN(PHI)**2*DM9~ 
X2TAYXRI*B12~ 
DX2ACYX14B12~ 
Yy20*(YYKI*B12#1/A12)- 
py2acyy1*B12- 
X3P*COS(PHI)*SIN(PHI)#S43~ 
DX3*COS(PHI)*SIN(PHT)*DM3~ 
YSTWSIN( PHI ) #28547 - 
DY3*SIN(PHI)#*24DM7+ 
‘TH3T*R34SINCPHI)##2*S41 1+ 
DTHSARSASIN(PHI)##24DM11 

“FSGNHAHIAWI4#2ABS (YIT) CONST 

-RMIMELW1##2#5 IN(WI#TS+PST1) 

~A22/A12#XM2#E24W 148245 IN(WI*TS+PSI2)) 
/CRMIAPSHHZ) -YIT. 


eee eee eee eee es 


RY1=( (YIT#DTY1)#24(DY1/P5)#82)4. 5 

YLI=24RYLASIN(PS*DTS/2) 

CALL PSIS(YIT, DY1, DTY1, PSY1,P5) 

CALL PHASE(Y1T, DY1,DTY1, PSY1,ANG, SGN1,SGN2,P5,DT5) 
yITC=¥1T 
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DYIC=DY1 

YIT=Y1T#SGNI*YLI*SIN(ANG) 

DY1=(DY1/PS+SGN2*YL1*COS(ANG) )*P5, 


c ITERATION 


6501 YBI=(¥iT#v1TC)/2 
DYB1=(DY14DY1C)/2 


CALL SIGN(DY1,SGN) 


DIYB1=(XB1*(COS(PHI)*SIN(PHI)*SM1+YXK1*B11)+ 
DXB1(COS( PHI) *SIN(PHI)*DM14CYX1#B11)+ 
YBI#(SIN(PHI)##2*SHS+YYK1*B1141/A11)+ 
DYBL*(SIN(PHI)##2#DMS+CYY14B11)+ 
‘THBL*RI*SIN( PHI) *#2#SN9+ 
DTHBI*RI*SIN( PHT) **2*DM9~ 
XB2AYXKI*B12~ 
DxB2*CYXI*B12~ 
‘YB2*(YYKI*B12+1/A12)~ 
DyB2*CYYI*B12~ 
XB3*COS(PHT)*SINPHI)*SM3- 
DXB3*COS( PHI) *SIN(PHI)*DM3~ 
YB3*SIN(PHI)*#2#S47 ~ 
DYB3*SIN(PHI)*#24DH7+ 
‘THBS*RS*SIN( PHI) *#24SH11+ 
DTHBS*R3*SIN(PHI)**2*DH11 

-#SGN*NO1#W 102 ABS (YB1)*CONST 

-KMIAELAW1**2*S IN(WIATS4PST1) 

~A22/A12#KM2*E2AW 1442S IN(WIATS+PSI2)) 
/ (SMTP 5442) -YB1 


eee eee eee eee eee ees 


RYBI=( (YITC+DTYB1)##2+(DY1C/P5)**2)**. 5 
YLB1=24RYBIASIN(P5*DTS/2) 
CALL PSIS(¥1TC, DYIC, DTYB1, PSYB1,P5) 
CALL PHASE(Y1TC, DY1C, DTYB1, PSYB1,ANG, SGN1, SGN2,P5,DT5) 
YBI1=YITC+SGNIAYLBI*SIN(ANG) 
DY1G/P54SGN2*YLB1*COS( ANG) )*P5 
‘ABS(Y1T) ~ABS(YB1) )*#2+(ABS(DY1) 
* _/P5-ABS(DYB1)/P5)**2 
TF(DLYL. GT. ACCR*(Y1T#*2+(DY1/P5)##2) )THEN 
YiT=B1 
pyl=pya1 
Go TO 6501 
ELSE 
GO TO 6502 
ENDIF 


6502 TS=TS4DTS 
6500 CONTINUE 
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DO 6600 LLL=1,1v2 


CALL SIGN(DY2, SGN) 


XI TAYXK1*B12~ 
* DxI*cYx14B12~ 

* YIT*(YYKI*B1241/A12)- 
* byI#CyYI#B12+ 

* X2TAYNKIMB22+ 

* Dx2acyxi#B22+ 

* YOT*(YYKI*B2241/A22)+ 
* 


Dy2*CYY14B22 
+#SGN*XM2#W142KABS (Y2T)*CONST. 
~XM2*EZAW 1828S IN(W1*T64PSI2) 
~ALL/AL24XM1*E1#W1#%24SIN(WI#T64PSI1) ) 

/(xMaeP6re2) -Y2T 


Crarepry2 )+42+(DY2/P6)m2)*. 5 
2*RY2*SIN(PO*DT6/2) 
CALL PSIS(Y2T, DY2, DTY2, PSY2,P6) 
CALL PHASE(Y2T, DY2,DTY2, PSY2, ANG, SGN1,SGN2,P6,DT6) 
yoT=¥2T 
Dy2c=py2 
YoT=Y2T+SGNI*YL2*SIN(ANG) 
DY2=(DY2/P6+SGN2*YL2*COS(ANG))#P6 


¢ ITERATION 


6601 yB2=(y2Te¥2TC)/2 
Dy24Dy2c)/2 


CALL SIGN(DYB2,SGN) 


= (-XBI*YXKI*B12~ 
* DXBI*CYX1*B12~ 

* YBI*(YYKINB1241/A12)- 
* pyBI*cyY1#B12+ 

* XB2*YXKI*B22¢ 

* DxB2*CYXI*B22+ 

* YB2*(YYKI®B2241/422)+ 
* 


byp2scyyien22 

-+SGN#XH2AW12#ABS (YB2)*CONST. 

12*EZAWLk24S IN(WIT6+PSI2) 

~ALL/AL2#XM1*E1*W19% 24S IN(W1*T64PSI1)) 
/(M2AP6HHD) -YB2 


RYB2=( (¥2TCHDTYB2 #24 (DY2C/P6)##2)%, 5 

YLB2=24RYB2#SIN(P6#DT6/2) 

CALL PSIS(¥Y2TC, DY2C, DTYB2, PSYB2, P6) 

CALL PHASE(Y2TC, DY2C, DTYB2, PSYB2, ANG, SGN1 ,SGN2, P6,DT6) 
YB2=¥2TC+SGNI*YLB2*SIN(ANG) 
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DY2C/P6+SGN2*YLB2*COS( ANG) )*P6 
‘ABS(Y2T) -ABS(YB2) )**2+(ABS(DY2) 
%* _/P6-ABS(DYB2)/P6)**2 
IF(DLY2. GT. ACCR#(Y2T¥2+(DY2/P6)*#2) )THEN 
yaT=¥B2 
DY2=DYB2 
GO TO 6601 
ELSE 
GO TO 6602 
ENDIF 


6602 ‘T6=T64DTS 
6600 CONTINUE 


DO 6700 LLLL=1,1Y3 


SM7=KO*(1-A*COS(WH*T7)) 
DM7=2T#2#( SM7#*XMS#XMI / (XM3+KH1) 9.5 


CALL SIGN(DY3,SGN) 
ory: 


~X1T*COS( PHT) *SIN(PHI)*SM1- 
DX1*COS(PHI)*SIN( PHI) *DM1- 
YIT*SIN(PHI)*#24S45~ 
DYISIN(PHI)**24D45~ 

‘THYT*RI*SIN(PHT) #24349 - 
DIHI*RI*SIN(PHI)*#2*DH9+ 
XST#(COS(PHI)#SIN( PHI )*SM34YXK2*B33)+ 
DX3#(COS( PHI)*SIN( PHI )*DM3+CYX2*B33)+ 
YST*(SIN(PHI)*#2#SH7-+YYK2*B334+1/A33)+ 
DY3*(SIN(PHI)*#2*DM74+CYY2#833) ~ 
‘THST*RI*SIN( PHI) **2*SMI1~ 
DTHS*RI*SIN(PHI)**2*DMI1- 

XG TAYRK2#D34~ 

DRGACYX2*B34~ 

YaT*(YYK2*B34+1/A94)~ 

DYGACYY2*R34 
+SGN#XM34W2*+#2*ABS(Y3T)*CONST 
-KSES*W2#%2#S IN(W2*T74PSI3) 

~AG4 /A34*KHGHEGHW2*H2*SIN(W2ATTHPSTS)) 
/OM3*P7##2)-Y3T 


eee eee ee ees 


YST4DTY3)*#24(DY3/P7) #82). 5 
RYSASIN(P7ADT7/2) 
CALL PSIS(¥3T, DY3, DTY3, PSY3,P7) 
CALL PHASE(Y3T, DY3, DTY3, PSY3, ANG, SGN1, SGN2,P7,DT7) 
Y3TC=¥3T 
DY3C=Dy3 
Y3T=¥3T#SGNIAYL3*SIN(ANG) 
DY3=(DY3/P74SGN2*YL3#COS(ANG) )*P7 
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c ITERATION 


6701 YB3=(¥3T#Y3TC)/2 
DYB3=(DY34DY3C) /2 


CALL SIGN(DYB3,SGN) 
prve: 


-XB1*COS(PHI)#*SIN(PHI)*SH1- 
DXB1*COS(PHT)#SIN(PHI)*DM1- 
‘YBI*SIN(PHI)**24SN5 ~ 
DYBI*SIN( PHT) #*2*DM5 ~ 
‘THBIARI*SIN(PHI)**2*SM9 ~ 
DTHBI*RI*SIN( PHT) *2*DMOF 
XB3*(COS( PHI) *SIN(PHI)*SM3+YXK2*B33)+ 
DXB34(COS( PHI) *SIN(PHI)*DM34CYX2*B33)+ 
YB3*(SIN(PHI)##2#SM7-+YYK2*B3341/A33)+ 
DYB34(SIN( PHI) *#2#DM74CYY2*B33)~ 
‘THB3*R3*SIN(PHI)**2*SM11- 
DITHB3*R3*SIN( PHI) **2*DM11- 
XBORYXK2AB34~ 
DXBOACYX2*B34- 

YBO*(YYK2*B34Q+1/A34)~ 
DyB4tCYY2*834 

-+SGN#XHS#W2#2#ABS (YB3)*CONST 

KHSAESAW2*#24SIN(W2#T7+PSI3) 

AGG / ABGAXMGAEGAWDHDRSIN(W2*T74PSI4)) 
/(XMS*P7##2) -YB3 


eee eee eee eee es 


(Y3TCHDTYBS)++2+(DY3C/P7)#*#2)*. 5 
A RYB3*SIN(PT#DT7/2) 
CALL PSIS(Y3TC, DY3C, DTYB3, PSYB3, P7) 

CALL PHASE(Y3TG, DY3C, DTYB3, PSYB3, ANG, SGN1, SGN2,P7,DT7) 
‘YB3=Y3TC+SGNI4YLB3*SIN(ANG) 
DYB3=(DY3C/P7+SGN2*YLB3*COS( ANG) )*P7 

DLY3=(ABS(Y3T) -ABS(YB3) )*#24(ABS(DY3) 

* /P7-ABS(DYB3)/P7)##2 

IF(DLY3. GT. ACCR#(Y3D%2+(DY3/P7)**2) )THEN 
YST=¥B3 
DY3=DYB3 
Go TO 6701 

ELSE 
GO 70 6702 

ENDIF 


67020 T7=T74DT7 
6700 CONTINUE 


DO 6800 LLLLI=1,1¥4 


CALL SIGN(DY4,SGN) 
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DTY4=( -X3TAYXK2*B34- 
DX3*CYX24B34~ 
Y3T*(YYK2*B3441/A34)- 
DYSACYYZ*B34+ 
XGTAYXKIARGG+ 
DXGHCYX2#BOA+ 
YOT#CYYK2#B4G41/864)+ 
DY4HCYYZ*BGG 
SGNHANGHW2#2*ABS (Y4T)*CONST 
~KHG*EG#W 244245 IN(W2*TB+PSI4) 
-A33/AIGAXMSAES*W2##24SIN(W2*TB+PSI3)) 
/CMGRPBERZ) -YOT. 


Pee e RE 


RYG=( (YOTHDTY4)#424(DY4/PB)i2) He, 5 
YLA=2*RY4#SIN(P8*DT8/2) 
CALL PSIS(Y4T, DY4, DTY4, PSY4,P8) 
CALL PHASE(Y4T, DY4 , DTY4,, PSY4, ANG, SGN1, SGN2,P8,DT8) 
yarTC=Yor 
DY4C=DY4 
YGT=Y4T+SGNI*YL4*SIN(ANG) 
DY4=(DY4/P8+SGN2*Y14*COS(ANG))*PB 


© ITERATION 


6801 YB4=(YaT+¥oTC)/2 
DYB4=(DY44DY4C) /2 


CALL SIGN(DYB4, SGN) 


DTYB4=(-XB3*YXK2*B34- 
DXxB3ACYX2*B34- 
YBS4(YYK2AB3G41/A34)~ 
DYBSACYY2*B34+ 
XBGAYNKZ#ROGE 
DXBGHCYX2*BG4+ 
YBO*(YYK2ABGG+I/AG4)+ 
DYBSACYY2*BG4 
FSGN#KMG*Y2A#2HABS (YB4)*CONST 
= KNGREGAW 2428S IN(W2*TB+PST4) 
A33/A3G#XM3*ES*W2*#2*SIN(W2*T84PSI3)) 
/CMGAPBHHD) -YBG 


errr ees 


RYB4=( (YOTCHDTYB4)#*24(DY4C/PB) HD) eH, 5 

YLBG=2*RYB4#SIN(P8*DT8/2) 

CALL PSIS(Y4TC,DY4C, DTYB4, PSYB4, P8) 

CALL PHASE(Y4TC, DY4C, DTYB4 , PSYB4, ANG, SGN1, SGN2, P8,DT8) 
‘YB4=Y4TC+SGNIAYLBG*SIN(ANG) 
DYB4=(DY4C/P8+SGN2*YLB4*COS( ANG) )*P6 

DLY4=(ABS(Y4T) -ABS(YB4) )*#2+(ABS(DY4) 

* —_/P8-ABS(DYB4)/P8)eA2 

TF(DLY4. GT. ACCR®(Y4T##2+(D¥4/P8)*#2) )THEN 
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YOT=¥B4 
DY4=DYB4 
Go TO 6801 
ELSE 
Go To 6802 
ENDIF 


6802 © T8=T8+DTS 
6800 CONTINUE 


Do 6900 NN=1, ITH 


‘SH9=KO*( 1-A*COS(WH*TI) 
DM9=21%2#( SHO*KMSAXMI / (XMS+K01) J, 5 


CALL SIGN(DTHI, SGN) 


DITHI=(X1T*RI*COS(PHI)**2*SM1+ 
DXI*RI*COS( PHI) ##2*DN 1+ 
YITARIASIN( PHI) ##24SHS+ 
DYIARI*SIN( PHI) *#2#DN5+ 
THUT#(RIS2#SHI*COS (PHI )*24KTL)+ 
DIHIARI*2*DM9ACOS( PHT) **2- 
‘THoreKTi- 

XOTARI*COS(PHI)*#2#SH3- 
DX3ARIACOS( PHT)*#2#DM3 - 
YST*RIASIN(PHI) #428847 ~ 
DYSARIASIN(PHI)##2*DH7+ 
‘THSTARIARS*SH11*COS (PHT) #2+ 
DTHS*RIARS*DM11#COS( PHI) #2 
+SGNAFRCI*2*WEI*D1/2) 
/(SV*P9*E2) -THIT, 


eee Tere Tees, 


(THIT#DTTH1)*#2+(DTH1/P9)**2)*. 5 
‘THLA=2#RTHI*SIN(P9*DT9/2) 
CALL PSIS(THIT, DTH1, DITH1, PSTH1, P9) 
CALL PHASE(THIT, DTH1, DTTHL, PSTHLL, ANG, SGN1, SGN2,P9,DT9) 
‘THITO=THIT 
DTHIC=DTHI 
‘THIT=THIT+SGN1*THLT*SIN(ANG) 
DTH1=(DTH1/P9+SGN2*THLI*COS(ANG))#P9 


ITERATION 


6901 THBI=(THIT#THITC) /2 
DTWB1=(DTHL+DTHIC) /2 


CALL SIGN(DTHBI, SGN) 


DITHBI=(XB1*RI*COS(PHI)*#2#SM1+ 
* DXBI*RI*COS(PHI)*#2*DM1+ 


149 
YBI*RIASIN( PHI) 2485+ 


DYBI*RI*SIN( PHT) **2*DMS+ 
‘THBI*(R1*2*SHI*COS (PHT) #24KT1) + 
DTHB1*R1*#2#DH9*COS(PHI)*#2- 
‘THB2*KTI- 
XB3*RIACOS(PHI)**2*SH3- 
DXB3*RI*COS(PHI)**2*DH3- 
YBSARINSIN(PHT )**24SM7~ 
DYBIARI*SIN( PHT) **2*DM7+ 
‘THBS*R1*R3*SH114COS( PHT) ##2+ 
DTHB3#R1#R3*DM11*COS (PHI) #2 

#SGNAPRCI*2#WE1*D1/2) 

/ (31962) -THBL 


eee es 


(THA TC+DTTHB1)*#24(DTHIC/P9) #2). 5 

2 RTHBI*SIN(PI*DT /2) 

CALL PSIS(THITC, DTHIC, DTTHB1, PSTHBI, P9) 

CALL PHASE(TH1TC,DTHIC, DTTHB1, PSTHBI , ANG, SGN1, SGN2,P9,DT9) 
‘THB1=TH1TC+SGN1*THLBI*SIN(ANG) 
DTHB1=(DTH1C/P9+SGN2*THLBI*COS(ANG))*P9 


DLTH1=(ABS(THIT) ~ABS (THB1) )*#2+(ABS(DTH1) 
‘* —_/P9-ABS(DTHB1)/P9)#2 
IF(DLTHL. GT. ACCR*CTHIT2+(DTH1/P9)**2) )THEN 
‘THIT=THB1 
Dri1=pTHB1 
Go TO 6901 
ELSE 
Go TO 6902 
ENDIF 


6902 T9=T94DT9 
6900 CONTINUE, 


DO 8000 NNN=1,ITH2 
CALL SIGN(DTH2, SGN) 


DrTH2=(-THIT*KT1+ 


* ‘THaT*KTI- 
* ‘TQLMSIN(WEX*T104X12) 
* 4SGNAFRC1*2#WE1*D1/2) 


/(32#P10*#2) -THAT 


RTH2=( (THAT#DTTH2 "2+ (DTH2/P10)*#2)#*, 5 
‘THL2=2*RTH2#SIN(P1O4DT10/2) 
CALL PSIS(TH2T, DTH2, DITH2, PSTH2, P10) 
CALL PHASE(TH2T, DTH2, DITH2, PSTH2, ANG, SGN1,SGN2, P10, DT10) 
‘Tu2Tc=TH2T 
DrH2c=pTH2 
‘TH2T=THAT+SGN1*THL2*SIN(ANG) 
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DTH2=(DTH2/P10+SGN2*THL2*COS(ANG))*P10 


c ITERATION 


8001 THB2=(THAT#TH2TC) /2 
DTHB2=(DTH2+DTH2C) /2 


CALL SIGN(DTHB2, SGN) 
DITHB2=(~THBI*KTL+ 
‘THB2*KT1- 
‘TQL*SIN(WEX*T104X12) 
+SGN*FRO1#*2*WE1*D1/2) 
/(324P10%*2)-THBZ 


seas 


RTHB2=(( (TH2TC+DITHB2) /100)*#2+(DTH2C/P10/100)*#2)**, 5*100 

‘THLB2=2*RTHB2*SIN(P10*DT10/2) 

CALL PSIS(TH2TC,DTH2C, DTTHB2, PSTHB2,P10) 

CALL PHASE(TH2TC, DTH2C, DTTHB2, PSTHB2, ANG, SGN1, SGN2,P10,DT10) 
‘THB2=TH2TC+SGN1*THLB2*SIN(ANG) 
DTWB2=(DTH2C/P10+SGN2*THLB2*COS(ANG))*P10 

DLTH2=(ABS(TH2T) ~ABS (THB2) )**#2+(ABS(DTH2) 

i /P10-ABS(DTHB2)/P10)*#2 

TRCDLTH2. GT. ACCR*(THETH*2+(DTH2/P10)*%2))THEN 
‘Tu2T=THB2 
DTH2=DTHB2 
Go To 8001 

ELSE 
Go 70 8002 

ENDIF 


8002 Tio=T10+DT10 
8000 CONTINUE 


DO 6100 NNNN=1, ITH 


sM1 
DML 


cOF(1-A¥COS(WH¥T11)) 
2 (SHI TRXMSAXMI /(XM34XMI) Yo. 5 


CALL SIGN(DTH3, SGN) 


DITHS=(X1TR3*COS (PHI )**2eSHI+ 
DXI*R3ACOS (PHI) **2*DM1+ 
YITARI*SIN(PHI) *24SH54+ 
DYI*R3*SIN(PHI)*#2*DN5+ 
THIT#RI®RS*SHO*COS (PHI) #24 
DTHI*RI*RI#DMI*COS(PHI)##2~ 
XST*R3ACOS( PHT) *2*SM3~ 
DX3*RI*COS( PHT) **2*DM3- 
YSTARS*SIN( PHI) "#2%SM7~ 
DY34R3¥SIN( PHI) **24DN7+ 


tee ee eee 
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‘THST*(R3**24SM1 1*COS( PHI) *#24KT2)+ 
DTHS*R3##2#DN1LACOS( PHI) ##2~ 
THGT*RT2 
“+SGN*FRO2*2*WE3*D2/2) 
/(33*P1182) “THAT 


teen 


RTH3=( (TH3T+DTTHS )**2+(DTHS/P11)*#2) 4, 5 
‘THLS=2*RTHI*SIN(PLI*DT11/2) 

CALL PSIS(THST,DTH3, DTTH3, PSTH3, P11) 

CALL PHASE(THST, DTH3, DITHS, PSTH3, ANG, SGN1, SGN2, P11,DT11) 


‘THSTC=THST 
DTH3C=DTH3 


‘THST=TH3T+SGNI*THL3*SIN(ANG) 
DTH3=(DTH3/P11+SGN2*THL3*COS(ANG) )*P11 


c ITERATION 


8101 THB: 
DTH: 


‘THSTHTHSTC) /2 
‘DTH3+DTH3C) /2 


CALL SIGN(DTHBS, SGN) 


DITHBS=(XB14R3*COS(PHI)*#2*SM1+ 
DXBI*R3*COS (PHI) ##2*DM1+ 
YBI*RS*SIN(PHI)*#2#SN5+ 
DYBI*R3*SIN(PHI)*#2#DN5+ 
‘THBIARI*RI*SH9*COS (PHI) #24 
DTHBIARI*RI*DM9*COS (PHI )*#2~ 
XB3*R3*COS(PHI)#2%SH3- 

DXB3*R3*COS (PHI )**2#DM3 - 
YB3*R3*SIN(PHI)#2*SM7~ 
DYB3*RIASIN(PHI)#*2*DM7+ 
‘THBS*(R34#2#SM11*COS(PHI)*#2+KT2)+ 
DTHBI*RI**2*DM11*COS(PHI)**2- 
‘THBG*RTZ 

#SGNAFRC2*2AWEI*D2/2) 

/(I3*P11**2) -THBS 


eee eee eee 


a 


(CTHSTC+DTTHB3)*#2+(DTH3C/P11)**2)*%. 5 
‘THLB3=24RTHBIASIN(PI1*DT11/2) 

CALL PSIS(TH3TC, DTIC, DITHBS, PSTHBS, P11) 

CALL PHASE(THSTC, DTH3C, DTTHB3, PSTHBS, ANG, SGN1,SGN2, P11,DT11) 
‘THB3=TH3TC+SGN1*THLBS*SIN(ANG) 
DTHB3=(DTH3C/P11+SGN2*THLBS*COS( ANG) )*P11 

DLTH3=(ABS(TH3T) -ABS(THB3) )##2+( ABS(DTH3) 

* /P11-ABS(DTHB3)/P11)##2 

TR(DLTH3. GT. ACCR*(THIT##2+(DTH3/P11)*#2) )THEN 
THST=THB3 
DTH3=DTHB3 
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Go To 8101 


ELSE 
Go TO 8102 
ENDIF 


8102 T11=T114DT11 
8100 CONTINUE 


DO 8200 NNNNN=1, ITH 
CALL SIGN(DTH4, SGN) 


DITHG=(~THST*KT2 
ATHGTARTZ 
#SGNAFRC2*2WE3*D2/2 
~TQ2*SIN(WEX*T124X14)) 
/ (3412642) -THOT, 


eee 


RTHO=( (THOT#DTTHS)#424(DTHG/P12)##2)%0, 5 

‘THLA=2*RTHG*SIN(P124DT12/2) 

CALL PSIS(THST, DTH4, DITHS, PSTHS, P12) 

CALL PHASE (THT, DIH4, DITH4, PSTH4, ANG, SGN1,SGN2,P12,DT12) 
‘THOTC=THOT 

4 

‘THAT=TH4T#+SGNI*THLA*SIN(ANG) 

DTH4=(DTH4/P124SGN2*THL4*COS(ANG) )*P12 


ITERATION 


8201 ‘THOTHTHOTC) /2 


DTH4+DTHOC) /2 


CALL SIGN(DTHB4, SGN) 


bras: 


~THB3*KT2 
“#THBGARTZ 
+SGNAPROZA2AWES*D2/2 
~TQ2*SIN(WEX*T124X14)) 
1 (aoe 122) ~THB4 


tee 


RTHB4=( (THOTC+DTTHBS #24 (DTHGC/P12)*#2)%*. 5 

‘THLB4=2*RTHBGASIN(P124DT12/2) 

CALL PSIS(TH4TC, DTC, DITHBS, PSTHBS, P12) 

CALL PHASE(THATC, DTH4C, DITHB4, PSTHBS, ANG, SGN1,SGN2,P12,DT12) 
‘THB4=TH4TC+SGNI*THLBG*SIN(ANG) 
DTHB4=(DTH4C/P12+SGN2*THLB4*COS(ANG))*P12 

DLTH4=(ABS (THAT) ~ABS(THB4 ) )**2+(ABS( DTH) 

‘* —/P12-ABS(DTHB4)/P12) #82 

TF(DLTHS. GT. ACCRA(THGT##2+(DTH4/P12)#2) ) THEN 
THOT=THBG 
DTH4=DTHBG 
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Go TO 8201 
ELSE 
Go TO 8202 
ENDIF 


8202 T12=T124DT12 
8200 CONTINUE 


6000 CONTINUE 


xA(Js41)=x1T. 

x2(Js41)=K2T 

XS(JH=KST 

X4(IS41| 

¥1 (341) 

¥2(J541)=¥2T 

Y3(JU#1)=¥3T. 

Y4(I541)=¥4T 

‘THA(JJ#1)=THIT: 

TH2(35+1)=TH2T 

TH3(J+1)=TH3T 

‘TMG (3341) =THET. 
DIS(JJ+#1)=RI*THI(JJ+1)4+RI*THS(JI+1) 
DGC JI+1)=(X1(IS+1)H2+Y (IIH) HEDIS 
DPCIS+1)=(X3( J+) 4243 (JIHL) HD) HH, 5 


7000 CONTINUE 


WRITE(30,919)N, GR, BD1, 11 ,WNO,WEX 
919 FORMAT(*RPM=", 116, 1X, "GEAR RAT. 


(= »1P4.1,1X,"BD=" 175.3, 
Lis", 1F6.3, 1X, 'WNO=! 


* 1x, y1F10.3, 1X, 'WEX=",1F10.3) 
WRITE(30,%*) 
WRITE(30,309) INCRE, ACCR 

309 FORMAT(2X,"INCR. ANGLE=",,1F11.8,5X,"AC. CR.=", 1F9.6) 
WRITE(30,*) 
WRITE(30, 319) 

319 FORMAT(4X, "TIME", 18X,"X1" , 15X, "SAMPLING? ) 
TIS=.0 


DO 311 JJ=1,JIS,IJIN 
WRITE(30,312) TJJ,X1(3I), JI 

312 FORMAT(IF11.7, 6X, 118.9, 10X, 116) 
‘TJJ=TIJ4DTIS*IIIN 

311 CONTINUE 


WRITE(30,#) 
WRITE(30,319) 

319 FORMAT(4X, ‘TIME , 18X,'X1", 15X, 'SAMPLING') 
TJJ=.0 
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pO 311 JJ=1,JJJ, JIN 
WRITE(30,312) TJJ,X1(JJ), JJ 

312 FORMAT(IF11.7, 6X, 1E18.9, 10X, 116) 
TII=TISHDTIIEIIIN 

311 CONTINUE 


WRITE(30,*) 
WRITE(30,*) 
WRITE(30, 329) 
329 FORMAT(4X, "TIME" , 18K, "X2",, 15X, "SAMPLING" ) 
TWJ=.0 
DO 322 JI=1,J33,IIIN 
WRITE(30,323) THJ,X2(JJ) JJ 
323 FORMAT(1F11.7,6X, 1E16.9, 10X, 116) 
TIJ=TISHDTIIFIIIN 
322 CONTINUE 


WRITE(30,#) 
WRITE(30,*) 
WRITE(30,339) 

339 FORMAT(4X, "TINE" ,16X,'X3", 15X, ‘SAMPLING" ) 
TH=.0 
DO 333 JJ=1,JJJ, JIN 
WRITE(30,334) T2J,X3(JJ), JJ 

334 FORMAT(IFL1. 7, 6X, 1E18.9, 10X, 116) 
‘TII=TISDTIICIIIN 

333 CONTINUE 


WRITE(30,*) 
WRITE(30,*) 
WRITE(30, 349) 

349 FORMAT(4X, "TIME" ,18X,'X4", 15X, "SAMPLING" ) 
TWJ=.0 
DO 344 JJ=1,JIJ,JJIN 
WRITE(30,345) TIJ,X4(JJ), JJ 

345 FORMAT(1F11.7,6X, 118.9, 10X, 116) 
TIJ=TIS4DTIIFIIIN 

344 CONTINUE 


WRITE(30,*) 
WRITE(30,#) 
WRITE(30, 359) 

359 FORMAT(4X, "TIME", 18X, "Y1", 15X, 'SAMPLING' ) 
TIJ=.0 
DO 355 JJ=1,JJJ,JJIN 
WRITE(30,356) TIJ,Y1(JJ),JJ 

356 FORMAT(1F11.7,6X, 1E18.9, 10X, 116) 
TIS-TISHDTIIEIIIN 

355 CONTINUE 


WRITE(30,*) 
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WRITE(30,*) 
WRITE(30,369) 
369 FORMAT(4X, "TIME", 18X, "Y2" , 15X, "SAMPLING" ) 
TIJ=.0 
DO 366 JJ=1,JIJ,JJIN 
WRITE(30, 367) TIJ,¥2(IJ), JJ 
367 FORMAT(1F11.7, 6X, 1E18.9, 10X, 116) 
‘TIJ=TISHDTIIEIIIN 
366 CONTINUE 


WRITE(30,*) 
WRITE(30,*) 
WRITE(30,379) 
379 FORMAT(4X, "TIME" , 18X,'Y3", 15X, ‘SAMPLING") 
TJJ=.0 
DO 377 JJ=1,IIJ,IIIN 
WRITE(30,378) TJJ,Y3(J),JS 
378 FORMAT(1F11.7,6X, 1E18.9,10X, 116) 
TIJS-TISDTIINIIIN 
377 CONTINUE 


WRITE(30, *) 
WRITE(30,*) 
WRITE(30, 389) 

389 FORMAT(4X, "TIME", 18X,"Y4", 15X, SAMPLING") 
TWJ=.0 
DO 388 JJ=1,JJJ, IJIN 
WRITE(30,387) TJJ,Y4(JI),IS 

387 FORMAT(1F11.7, 6X, 1E18.9, 10X, 116) 
TIJ=TISDTII*IIIN 

388 CONTINUE 


WRITE 30, #) 
WRITE(30,*) 
WRITE(30,397) 

397 FORMAT(4X, "TIME", 18X, "THI", 14X, ' SAMPLING") 
‘TH=.0 
DO 399 JJ=1, JIS, IJIN 
WRITE(30,398) TIJ,THI(II) JI 

398 FORMAT( 1F11.7, 6X, 1E18.9, 10X, 116) 
TIJ-TIIHDTIINIIIN 

399 CONTINUE 


WRITE(30,*) 
WRITE(30,%*) 
WRITE(30, 409) 
409 FORMAT(4X, ‘TIME, 18X, 'TH2", 14X, "SAMPLING" ) 
TH=.0 
DO 400 JJ=1,IJ3,IIIN 
WRITE(30,401) TJJ,TH2CII) JT 
401 FORMAT(1F11.7,6X, 1618.9, 10X,116) 
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‘TIJ=TII4DTIIAISIN 
400 CONTINUE 


WRITE(30,*) 
WRITE(30,*) 
WRITE(30,419) 

419 FORMAT(4X, ‘TIME, 18X, "TH3" ,14X, "SAMPLING" ) 
Ti. 


DO 411 JJ=1,JJJ, IJIN 
WRITE(30,412) TIJ,TH3(II) II 
412 FORMAT(1F11. 7, 6X, 1618.9, 10X, 116) 

TJJ=TISDTSIEIIIN 
411 CONTINUE 


WRITE(30,*) 
WRITE(30,*) 
WRITE(30, 429) 

429 FORMAT(4X, "TIME" , 18X, "THA" , 14X, " SAMPLING' ) 
TIJ=.0 
DO 422 JJ=1, JJ, IJIN 
WRITE(30,423) TJJ,TH4(IJ) JF 

423 FORMAT( 1F11. 7, 6X, 1618.9, 10X, 116) 
‘TIS=-TIIADTII*IIIN 

422 CONTINUE 


WRITE(30,*) 
WRITE(30,719) 
719 FORMAT(4X, 'TIME" , 18X, "DP", 15X, "SAMPLING" ) 
TIJ=.0 
DO 711 JJ=1,J4J,IIIN 
WRITE(30, 712) TIJ,DP(II) JT 
712 FORMAT(1F11.7, 6X, 1E18.9, 10X, 116) 


TISSTIS#DTIIASIIN 
711 CONTINUE 
WRITE(30,*) 
WRITE(30, 819) 
819 FORMAT(4X, "TIME" , 18x, "DG" , 15x, SAMPLING’ ) 
TyJ=.0 
DO 811 33=1,JJJ,JJIN 


WRITE(30,812) TJJ,DG(JJ), JJ 

812 FORMAT(1F11.7, 6X, 1£18.9, 10x, 116) 
TIJ=TIDTIIAIIIN 

611 CONTINUE 


WRITE(30,*) 
WRITE(30,918) 

918 FORMAT(4X, "TIME" , 18X, "DTS" , 15%, "SAMPLING ) 
TWJ=.0 
bo 911 JJ=1,5IJ, IJIN 
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WRITE(30,912) TJJ,DTS(JJ), JI 
912 FORMAT(1F11. 7, 6X, 1E18.9, 10X, 116) 
TIJ=TIJ4DTIICIIIN 
911 CONTINUE 


STOP 
END 


e990 


EIGENVALUES * 
* 
CALCULATES EIGENVALUES FOR TWELVE VARIABLES * 
* 


IMPLICIT REAL*® (A-H) 
IMPLICIT REAL*8 (0-2) 


REAL 11, 12,J1,J2,J3,J4,NU, 121, 122,KT1,KT2, INCRE, JEFF 
REAL L1,1.2,13,14,15,16,K0,K, KOAVE 


COMMON WE1,WE2,WE3,WE4,D1,D2,L1,L2,14,15,GR 
COMMON BL1,BL2,BD1,BD2,C,W1,W2,N 


CHARACTER*B FN30,FT30 
FN3O="EIGEN' 
FT30="DATA’ 

CALL cMs 


*C'FILEDEF 30 DISK '//FN30//' '//FT30//' * ( LRECL 255 RECFM V" 


GIVEN PARAMETERS-. 
‘D'S FOR SHAFTS, R'S FOR ROTORS 


CONV=180/3. 14159 


158 


wD 


E=30000000. 
NU=.292 

G=E/(2*(14NU)) 
‘TQ=WI*1000. 

N=47750 

WI=N#3. 14159%2/60 
W2=GREWI 

WE1=3. 14159*RI4H2WI*, 282 
WE2=WEL 

WE3=3. 141594RS#H2AWI*, 282 
WEO=WES 

XM1=WE1/CONVL 
XM2=WE2/CONV1 
XM3=WE3/CONVL 
XMG=WE4/CONVI 


14 1S9*DIHHG/64, 
M61 59#D2H4/64 


“INFLUENCE COEFFICIENT: 


ALI=1243/(6*E*11) 
A12=(LU*L2662) /(GHEXT1) 
AQQ=L1*#2#(L 142412) /(3*E*I1) 
AS3=15%43/(6*E#12) 
A34=(LAKLS #42) / (G*EXT2) 
AGG=LA th Dk( Lit 2*LS) / (3*EHT2) 


BL1=(LU#24(L14L2*2)#42) /(L14L2)**D 
B12=(2AL2HH2) /(L14L2) #82 

B22=B12 
B33=(L4th2+( LAL #2 )H2) /(LA4L5 42, 
B34=(2AL5##2) /(I44L5 #2 

BA4=B34 


121=3, 14159*D14#4/32. 
122=3.14159*D2*44/32 
KTISIZ18G/(L14L2) 
KT2=1224G/(LA4L5 )*GRIHZ 
KO=WI#1.2*1000000 
K=KOARIM#24COS(PHI)**2 


TARIMH2/2. 
1O*RIAHD/2 

SARIN /ZAGREMD 
IMGARGIAD/2AGRIRD 
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160 
WEX=WNO 


2=20 
wMezewL 


CALL SUBROUTINE FOR K'S & C'S DUE TO BEARING- 


CALL BRNG(S1, $2, XXK1,XXK2, XYK1,XYK2, YNK1, YXK2, YYK1, YYK2, 
* CXX1, CXX2, CXY1, CXY2, CYN1, CYX2,CYY1,CY¥2) 


(COS( PHI) **2*KO+XXKI*B1141/A11) /XML)**. 5, 
(CAXK1*B2241/A22) /XM2) 5 

(COS ( PHI)**2*KOFXXK2*B33+1 /A33) /XM3)#, 5 
CORK 2*BOGET/AGG) /XMG) A 5 
(SIN(PHI)**2*KO+YYRI*B1141/A11) /XML)##.5 
CYYK1*B2241/422) /XM1 #5 
(SIN(PHI)*#2*KO+YYK24B3341/A33) /XM3)*.5 
CYYR2*B4G+1/A44) /XMG)#%. 5 
(RTLAR1**2#KORCOS(PHT)*#2) / JL). 5 
WNaO=(KT1/J2)##.5 

WNI1=( (KT2+R3*#2#KORCOS( PHI)**2) /J3)*. 5 
WN12=(KT2/J4)*. 5 


ESSSSSS555 


SWN1=( (COS( PHI) 2*KO+XXK1*B1141/A11) /XML)**. 5 
SWN3=( (COS ( PHI) *#2*ROFXXK2*B33+1 /A33)/XM3)*. 5 
SWN5=( (SIN( PHI)*#2*KO+YYKI*BII+1/A11)/XM1) "8.5 
'SWN7=(( SIN( PHI) *#2*KO+YYR2*B3341/A33) /XM3)*. 5 
SWN9=(RI#24KORCOS (PHT) **2/J1)**. 5 
R3®R2AKORCOS (PHI )##2/13)4.5 


© UNCOUPLED NATURAL FREQUENCIES 


UNF1=(KT1/(J1452/(14I2)))**. 5 
INF1/1.414213, 

KT 2/(J3*34/(I34I4)) #5 

“UNF2/(2##.5) 


KOMCOS(PHI )*#2/ (XMI#XM3 /(XMI4KM3))):4.5 
KO*SIN(PHI)*82/(XMIAXM3/(XMI4KH3) ) 24.5 
(UNF35**2+UNF37#%2)/2)#*, 5 
UNF4=(1/(XHIAIL) 5 

UNFS=( (48*E#12)/((L5+L6)#3) /XM3)**. 5 
1/(x2#A22) Joe 5 

1/(XMI*A12) 5 

UNE7=( (6*EA12)/(2AL4G#H 24 (LO4L5+1.6)) /XMS) 5 
UNFB8=((34E*11) /(L1¥#3) /XM2) 5 
UNF9=((3*E*12) / (L443) /XMG)8 5 


161 


CONF 2#XM1FUNFS##2#XH3) /(XMI4KNS) J. 5 
UNF 1=(30XK1B12 /XM1)40, 5 
UNF111=(XXK1#B22/KH2), 5 
UNF12=(YYKI*B11/XM1) #5 
UNF121=(¥¥KI*B22/XM2)#+, 5 
UNF13=(XXK2#B11/43) #5 
UNF131=(10KK2*B22/XM4) #5 
UNF14=(YYK2#B11/X012)*.5 
UNF141=(YYK2#B22/XM4) #5 

UNF1S=( (UNFLI®#24UNF12#42)/2)9#.5 
UNF1S1=( (UNFLL1**24UNF121%%2) /2), 5 
UNF16=((UNF13**24UNF14%42)/2)44.5 
UNF161=((UNF131%#24UNF141##2) /2)#*. 5 


SUNF1=(2*34E411/L2403/XM1)é, 5 


WRITE(30,11)8D1,BE1,11,12,L3,R1 

211 FORHAT(' BDI" 1F7.4,2X, "BEE s1F7.4,2X,'L1=",1F7.4,2X, 

1F7.6,2X,"L3=" ,1F7.4,2K, "R12", 1F7-4) 
“arte, 112) GR, XM1/iM2 ,KO, K 

112. FORMATCAX, ' 1 1F5.2,2K, MASS RATIO=" ,1F4. 2 
* 1X, "KO=" 1611.6, 1X, 'K> 
WRITE(30, 113)S1,82,KT1,KT2 

113 FORMAT(IX, 'S1=", 1F10.4, 2X, 'S2=", 1F10.4,2X,'KTI=",1E14.6,2X,"KT2=", 
* 1E16..6) 
WRITE( 30, 114)XXK1, XXK2, XYK1, XYK2 

114 FORMAT( "XK 
* *XYK2=" , 1E11.4) 
WRITE(30, 15)YKKL, YXK2, YYKI, Y¥K2 

15 FORMAT( 'YXK1=" 1611.4, 2X, 'YXK2=" ,1E11.4, 2X, 'YYKIS", 1E11.4,2X, 
* "YyK2=", 1811.4) 
WRITE(30, 116)CXX1,,CXX2, CXY1, CXY2. 

116 FORMAT(' CXX1=" ,1B11.4,2X, ' CXX: 
* *Cxy2=", 1611.4) 
WRITE(30, 117)¢vx01, C¥x2, cy¥1,cv¥2 

TELL.4 2X, "CYX: 

»1B11,4) 
WRITE(30, 119)B11, 812,822 

119 FORMAT(1X, 'B11=" , 1F10.3, 3X, 'B12=", 1F10,3, 3X, 'B22=" ,1F10.3) 
WRITE(30,919)1/A11, 1/412, 1/422 

919 FORMATCIX, "1/A11=",1F15.5,3X,"1/A12=" , 1F15.5,3X, "1/42; 
WRITE(30, 120)XM1, xH2, XM 4 

120 FORMAT(1X, 'Mi=", 1F8.3,3X,'N2=" , 1F8.3,3X, 'M3=", 1F6.3, 3K, "MG: 
* ) 
WRITE(30, 190).1,J2,J3,54 

190 FORMAT(IX, 'J1="> 1F7.3, 3X, "J: 
* ) 
WRITE(30,*) 
WRITE(30, 123)WN1,WN2,WN3,WNG 

123 FORMAT(AX, 'WNI=" ,1FL1.3, 3X, "WN2=", 1F11.3,3X, 'WN3=", 1F11.3, 3%, 
* 1Fi1.3) 


»1ELL.4, 2X, "XYKI=", 1611.4,2X, 


»1E11.4,2X, "CXYI=",1E11.4,2X, 


*1E1L.4, 2X, 'CYYI=", 1611.4,2X, 


y1F15.5) 


s1F8.3 


1F7.3,3X, "335" ,1F7.3,3X,"J4=", 187.3 


ao 


eeee 
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WRITE( 30, 124)WNS,WN6,WN7,WNS. 


y1FL1,3,3X,'WN7=",1F11.3,3X, 


WRITE(30, 175)WN9,WN1O,WN11,WN12 
175 FORMAT(IX, 'WN9="  1F11.3, 3X, "WN1O: 
* *wN12=", 1F11.3) 
WRITE(30,*) 
WRITE( 30, 129)WEX,WNO, WM, JEFF 
129 FORMAT(1X, 'WEX=" 1F10.3,2X, 'WNO=" , 1F10.3,2X, "WH=", 1F10.3,2X, 
* *JEFF=" , 1F8.5) 
WRITE(30, 139)WEX/WNO, WEX/WM, WEX/ (21) ,W1,W2 
139 FORMAT( 'WEX/WNO=" , 1F4. 2, 1X, 'WEX/WH=" , 1F7.3, 1X, 'WEX/(Z#W1)=" 


1F11.3,3X, 'WNLI=" ,1F11.3,3X, 


* »1F7.3, 1X, 'WE1s", 1FB8.2, 1X, 'WE2=", 1F8.2) 
WRITE(30,*) 
WRITE(30, 501)UNFL 

501 FORMATCIX,"UNF1 (SHAFT TORSION) =, 1F11.3) 
WRITE( 30, 551)UNFIE. 

551 FORMAT(1X, 'UNFLE(EQV.TOR.(NEG.GEAR INER,1)=",1F11.3) 
WRITE(30, 502)UNF2 

902 FORMAT(IX, 'UNF2 ( 2ND SHAFT TORSION) =, 1F11.3) 
WRITE(30, 592)UNF21 

592 FORMAT(IX, 'UNF21(EQV. TOR. (NEG.GEAR INER,2)=",1F11.3) 
WRITE(30, 503 )UNF3 

503 FORMAT(1X, 'UNF3 (TORSIONAL MESH) =", 1F11.3) 
WRITE(30,513)UNF35 

513 FORMATCIX, 'UNF35 (LINEAR MESH, x) =", 1F11.3) 
WRITE(30,523)UNF37 

523 FORMATCIX, 'UNF37 (LINEAR MESH, Y) =',1F11.3) 
WRITE(30,533)UNF39 

533 FORMATCIX, "UNF39 (COMBINED LINEAR MESH)  =',1F11.3) 
WRITE(30, 504)UNF4 

504 FORMATCIX,'UNF4 (BENDING. SIMPL SUPP.) =",1F11.3) 
WRITE(30, 505 )UNFS 

505 FORMAT(1X,'UNFS (BENDING. SIMPL SUPP.) »1F11.3) 
WRITE(30, 506 )UNF6 

506 FORMAT(1X, 'UNF6 (BNDNG. SMPL SUP OVERHNG.)=",1F11.3) 
WRITE(30,526)UNF69 

526 FORMAT(1X, "UNF69(BNDG. CTF SMPL SP OVHNG. )=",1F11.3) 
WRITE(30,507)UNF7 

507 FORMATC1X, 'UNF7 (BNDNG. SMPL SUP OVERHNG. )=",1F11.3) 
WRITE(30,508)UNF8 

508 FORMATCIX, "UNF8 (BNDNG CANTILEVER) y1F11.3) 
WRITE(30, 509 )UNF9 

509 FORMAT(1X, 'UNF9 (BNDNG CANTILEVER) »1F11.3) 
WRITE(30, 510)UNF10 

510 FORMAT(1X, 'UNF10 (BNDNG, COMBINED) »1F11.3) 
WRITE(30,511)UNF11. 

511 FORMAT( 1X, 'UNF11 (BEARING SUPPORT, X1,1ST)=",1F11.3) 
WRITE(30,611)UNF111 

611 FORMAT(1X, 'UNF111(BEARING SUPPORT, X2,1ST)=",1F11.3) 


eecaneaa 


512 


612 


513 


613 


514 


614 


515 


615 


516 


616 


901 


902 


903 


905 


904 


WRITE(30,512)UNF12 
FORMAT(1X, 'UNF12 (BEARING SUPPORT, 
WRITE(30,612)UNF121 

FORMAT( 1X, 'UNF121(BEARING SUPPORT, 
WRITE(30,513)UNF13, 

FORMAT(1X, 'UNF13 (BEARING SUPPORT, 
WRITE(30,613)UNF131 


4F11.3) 


FORMAT( 1X, 'UNF131(BEARING SUPPORT, 
WRITE(30,514)UNF14 

FORMAT(1X, 'UNF14 (BEARING SUPPORT, Y1,2ND)=",1F11.3) 
WRITE(30,614)UNF141 

FORMAT( 1X, 'UNF141(BEARING SUPPORT, ¥2,2ND)=",1F11.3) 
WRITE(30,515)UNF15, 

FORMAT(1X, 'UNF15 (BRNG SUPPORT, 1, 1CMBINED)=" , 1F11.3) 
WRITE(30,615)UNF151 

FORMAT( 1X, 'UNF151(BRNG SUPPORT, 1, 2CMBINED)=", 1F11.3) 
WRITE(30,516)UNF16 

FORMAT(1X, 'UNF16 (BRNG SUPPORT, 2, 1CMBINED)=" , 1F11.3) 
WRITE(30,616)UNF161 

FORMAT( 1X, 'UNF161(BRNG SUPPORT, 2,2CMBINED)=",1F11.3) 


WRITE(30,*) 
WRITE(30,901) 


FORMAT(1X, "4" SIMPLIFIED SYSTEM *##") 

WRITE(30,*) 

WRITE(30, 902)SWN1, SWN3, SHNS 

FORMAT( 2X, "SWNI=", 1F11.3,3X, "SWN3=", 1F11.3,3X, 'SWNS=" 


WRITE(30, 903)SWN7, SWN9, SWNI1 
FORMAT(2X, "SWN7=" , 1F11.3, 9X, 'SWN9=", 1F11.3,3X, 'SWNIJ: 


WRITE(30,*) 

WRITE(30,905)UNF3 

FORMATCIX, "UNP3. (TORSTONAL MESH) 
WRITE (30, 904)SUNF1 

FORMAT( 1X, "SUNF1 (DOUBLE CANTILEVR, x) 


»1F11.3) 


1F11.3) 


STOP 
END 
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»AFLL.3,3X) 


*1F11.3,3X) 


eanaee 


eee chic 


eee 
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ee 


‘* THIS SUBROUTINE IS TO CALCULATE THE BEARING RELATED PARAMETERS * 


* SUCH AS SOMMERFELD NUMBER, STIFFNESS AND DAMPING COEFFICIENTS * 


* 


{eccentrics 


SUBROUTINE BRNG(S1,S2,XXK1 , XXK2, XYK1, XYK2, YXK1, YXK2, YYK1, YYK2, 
* CXX1, CXX2, CXY1,CXY2, CYX1, CYX2,CYY1,C¥¥2) 


IMPLICIT REAL*8 (A-H) 
IMPLICIT REAL*8 (0-2) 


REALM8 S(2, 500), SN(2),RS(2),RJ(2), BL(2), BD(2) 
REALMS PA(2),DTEM(2, 500), TMA(2, 500) ,VISCA(2,500) 
REAL*8 PS(2, 500), PVISCA(2,500) 

REAL#8 XXKP(2), XYKP(2) , YXKP(2), YYKP(2),TOP(2) 
REAL*8 CXXP(2), CXYP(2),CYXP(2) ,CYYP(2),R(2) 


REAL 11,12,13,14,15,16 
COMMON WE1,WE2,WES,WE4,D1,D2,11,L2,14,L5,GR 
COMMON BL1,BL2,BD1,8D2,C,W1,W2,N 


OIL IS SAE20----- 


WEL4WE243, 14159/48D1*428(L142412)*, 282)/2 
WEI4WEGHS. 14159 /4*D2**24( 1442415 J, 282)/2 
RS(1)=N/60. 
RS(2)=GRARS(1) 
RJ(1)=01/2 
RI(2)=2/2 
BL(1)=BL1 
BL(2)=BL2 
D(1)=BD1 
Bo(2)=Bb2 
‘TMI=110. 
PAC1)=WRT1/(D1*BL1) 
PAC2)=WRT2/(D2*BL2) 
VISCO=1. 36/1088 
PVISCO=VISCO*100000 


‘TOP(1)=. 8003+. 27834R(1) 
‘TOP(2)=. 8003+. 2783#R(2) 


‘DO LOOP FOR SOMMERFELD NO- 


bo 115 


DTEM(I, 1)=.0 
DO 111 J=1,499 


‘TMACL, J)=THL4DTEM(T,J)/2 

VISCACT, J)=VISCO*EXP(1271.5/(THA(I,J)+95)) 
PVISCA( I, J)=PVISCO*EXP( 1271.5/(THACT, J}+95)) 

S(1,J)=VISCA(T, J)*RS(I) /PACT)*(RI()/C)2 
PS(,J)=PVISCA(T, J)*RS( I) /PA(1)*(RI(1)/C)**2 


TF(R(L). GE. .5)THEN 
IF(PS(1,J). GT. 15000)THEN 


DTEN(1, J#1)=84. 2989/(R(1)**, 08167)#PS(1,J)**(. 8554+ 
* + 08787*R(I))*PACT) /JNC/100000%#TOP(I) 
ELSE 
DTEM(1,J+1)=42..0097/(R(I)**.4146)#PS( 1, J)**(. 6869- 
* = 16*R(1) )*PA(T)/JINC/100000#T0P(T) 
ENDIF 
ELSE 


IF(PS(I,J). GT. 15000)THEN 
DIEM( 1, J+1)=69..4842/(R(1)**. 3075)*PS(I,J)*#( 8003+ 


* -2783*R(1) )*PA(1) /JNC/100000%*T0P( 1) 
ELSE 
DTEM(T, J+1)=38. 7604/(R(T)#*.5307)*PS(I,J)*#(.7322- 
* + 2505*R(1) )*PACT) /JNC/100000*#T0P(1) 
ENDIF 
ENDIF 


‘TMACL, J41)=THT4DTEM( I, J#1)/2 

VISCAC I, J+1)=VISCO*EXP( 1271.5/(THACI, J#1)495)) 
PVISCA( 1, J+1)=PVISCO*EXP( 1271.5/(TMA(T, J+1)495)) 

S(1,J#1)=VISCA(T, J41)#RS(T) /PACT)#*(RI(L) /0)##2 
PS(I,J#1)=PVISCA(T, J+1)*RS( 1)/PA(I)#(RJ(1)/C)*#2 


TF(ABS(PS(T,J#1)-PS(I,J)). LT. 2000. )THEN 
SN(1)=PS(1,J+1)/100000. 
Go TO 115 

ELSE 
Go To 111 

ENDIF 


111 CONTINUE 
115 CONTINUE, 


DO LOOP FOR K's & C'S 


165 


166 
DO 222 K: 


Be 
IF(R(K). GE. .5)THEN 


IF(SN(K). GE. 1.0)THEN 
XXKP(K)=2. 3258-1, 2124R(K)-.3413*SN(K)+, 3436#SN(K)**R(K) 
XYKP(R)=8. 1515*R(K) #1. 14424SN(K)*#(. 3874.47 174R(K)) 

2. 2356-12. 4891*R(K)+1. 2669*SN(K)+. O756*SN(K) #2. 

#2. 2224*SN(K)*R(K)-10.9395#SN(K)"#R(K) 

2178*R(K)*4(~. 1386) *SN(K)##( ~. 0376+. 0383#R(K)) 


IF(SN(K).LT. .15)THEN 
XXKP(K)=.5979*R(K)##( 1.0181 )*SN(K)#( ~. 88634. 19274R(K)) 
XYKP(K)=2. 0501#R(K)*4( ~. 2127) #SN(K)*#( ~_ 37134. 1476#R(K)) 
YERP(K)=- 481641. 7006#R(K) ~.9335#SN(K)411. 694SN(K)#42 

* 716. 3368*SN(K)*R(K)42.2198*SN(K)*AR(K) 
YYKP(K)=2.5181#R(K)*#( ~. 3236)*SN(K)**(.4954-. 4OO7#R(K)) 

ELSE 
XXKP(K)=1. 125 1R(K)**( ~..6746)*SN(K)*#(~, 81794. 46914R(K)) 
XYKP(K)=7.5105*R(K)*#. 9778#SN(K)#H(~. 558441 .01314R(K)) 
YRRP(K)=1.46342..044R(K) 1. 29#SN(K)41.053#SN(K) #2 

* 712, 272#SN(K)*R(K) 44. 378*SN(K)*#R(K) 
YYKP(K)=2,22024R(K)*#(~. 1794 )#SN(K)**(3145-.27714R(K)) 

ENDIF 
ENDIF 


IF(SN(K). GT. 1.0)THEN 
CRXP(K)=26. 424-25 .646#R(K)~11.235#SN(K)~.04G*SN(K)#H2 


* 443, 39 1*SN(K)*R(K)~16. 9474SN(K)##R(K), 
CXYP(K)=1. 35341. 109#R(K)+. 3764SN(K)+. 027#SN(K)42 
* ~1. 341*SN(K)*R(K)+, 510*SN(K)R(K) 
GYXP(K)=4..059-2.0154R(K) ~.231*SN(K)-. OO4*SN(K)**2 
* +1. 839*SN(K)*R(K)=1, 450#SN(K)H#R(K) 
CYYP(K)=4. 735~5.019#R(K)~7. 506#SN(K) -. OB9#SN(K) #42 
* #26. 0B9*SN(K)*R(K) -2.897#SN(K)##R(K) 


ELSE 
TF(SN(K). GT. .25)THEN 

- 357-2. 298*R(K)~17, 2948SN(K)43. 521SN(K) #2 

* +19. B92SN(K)*R(K)45. 635#SN(K)**R(K) 

464141. 82048R(K) ~. 6852#SN(K)+. 3214#SN(K) #2 


* 2. 1612#SN(K)#R(K)+2, 2092#SN(K)**R(K) 
461942. 02404R(K) =~. 5145*SN(K)+,4168#SN(K)##2 

* ~2. 1685*SN(K)*R(K)+1. 98 17*SN(K)*#R(K) 
7.852410.416*R(K)~15. 1858SN(K)43, 567#SN(K)*2, 

* +10. 084#SN(K)*R(K)+13. 869#SN(K)**R(K) 


ELSE 
IF(SN(K).GT..05)THEN 
CXXP(K)=26..402~19. 317*R(K)~5. 201*SN(K)+55. 9454SN(K)##2 

* #17. 1024SN(K)*R(K) 29. 361*SN(K)*#R(K) 

CRYP(K)=1, 7126+. 1870#R(K)-9. 7244#SN(R)410. O934KSN(K)*RD 

* +5. 0982*SN(K)*R(K)42. 8168*SN(K)i"#R(K) 
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CYXP(K)=1.6125~. 201*R(K)-5, 4768*SN(K)+. 1344#SN(K)#82 
* +46. 7569*SN(K)*R(K)+2. 2069%SN(K)*R(K) 

CYYP(K)=1,587-1.083*R(K)~18. 793*SN(K) +13. 865*SN(K)*#2 
* $30. 244SN(K)#R(K)+2. 02 4SN(K)MER(K) 

ELSE 

OXXP(K)=-8.09+30. 59#R(K)~1647, 56#SN(K)47078. 46*SN(K)##2 
* +4739. BO*SN(K)*R(K)4219. 81SN(K)"R(K) 

CXYP(K)=~21. 11427. 88*R(K)~705 . 62*SN(K)42327. 284SN(K) #2 
* +353. O2*SN(K)*R(K)4133.45*SN(K)##R(K) 

CYXP(K)=-22. 24428. 43*R(K) ~672*SN(K)42129.. 06*SN(K)##2 
* 4337. 77#SN(K)*R(K)+132. 83*SN(K)*#R(K) 

CYYP(K)=-2.87245. 397#R(K) ~204. 892*SN(K)4833, 2230SN(K)##2 
* #106. O1*SN(K)AR(K)431. 329#SN(K)*#R(K) 

ENDIF 

ENDIF 

ENDIF 
ELSE 


TE(SN(K). GT. 1.0)THEN 
+ 27024R(K)*#( ~ . 5053) *SN(K)**(~. 8199+. 97234R(K)) 
-TA7TAR(K) +. 64O1¥SN(K)**( 2659+. 8173#R(K)) 
YXKP(K)=27 92-32. 71#R(K)=1. 21*SN(K)=.404#SN(K)**2 
* #19. 36*SN(K)*R(K) 22. 18SN(K)R(K) 
YYKP(K)=3. 1903*R(K)*, 3860*SN(K)#*(. 1813-. 3994#R(K)) 
ELSE 
IF(SN(K).LE. 15) THEN 


XXKP(K)=. 35324R(K)*#(~1.7179) *SN(K)**(~, 2589-1.09224R(K)) 
XYKP(K)=1.91654R(K)*#( ~ .2674)#*SN(K)#*( ~, 6257+. 6357#R(K)) 
YXKP(K)=16. 74-26. 534R(K)+13. 63*SN(K)~38. 7O*SN(K)#2 
* #43, 72*SNCK)*R(K) ~20. 12*SN(K)*#R(K) 
YYKP(K)=7 455 14R(K)##( 1, 1975) *SN(K)*#(~..611741. 8354*R(K)) 
ELSE 
XXKP(K)=1.1897#R(K)##( 
XYKP(K)=3. 81028R(K)##( 
IF(SN(K). LE. .5)THEN 
YXKP(K)=4.4333~6. 0498*R(K)+3, 9980*SN(K)-4. 21338SN(K)**2 
+. 8B4B*SNCK)AR(K) ~6. 692 1*SN(K)R(K) 


6076 )*SN(K)*#(~.7152+. 3053*R(K)) 
02386) *SN(K)##(~.44757+, 85571#R(K) ) 


ELSE 
YXKP(K)=-7..027410. 473#R(K) ~1. 964SN(K)42.039%SN(K)#2 
* ~194675*SN(K)*R(K)+8, 95*SN(K)*AR(K) 
ENDIF 


YYKP(K)=3.417 14R(K)**. 45O7#SN(K)*#(~, 2406+, 8095##R(K)) 
ENDIF 
ENDIF 


IF(SN(K). GT. 1.0)THEN 
CXXP(K)=33.5977-36. 630#R(K)~6. 205#SN(K)-. 915*SN(K)*H2 
#38. 34G#SN(K)*R(K)=20. 881*SN(K)AAR(K) 
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CRYP(K)=4, 872-5. 6578R(K)-. 391SN(K)+. 007#SN(K)##2 


* ++.576#SN(K)*R(K)+. 2018SN(K)*AR(K), 
CYRP(K)=5. 145~4.047#R(K) ~. 577#SN(K) ~. O16*SN(K)**2 
* #2. T4dSN(K)#R(K) “1. 615*SN(K)MER(K) 
CYYP(K)=3..657-7..447#R(K)-3.246%SN(K)+. 1744SN(K)**2 
* 411. 9324SN(K)*R(K)#1. 95 1SN(K)"ER(K) 


ELSE 
TE(SN(K).GT. .10)THEN 
GRXP(K)=45. 534-57. 339#R(K)+5. 652*SN(K)-7. 791*SN(K)**2 


* #61. O6*SN(K)*R(K)-38 .02#SN(K)"#R(K) 
CXYP(K)=8. 035-10. 785*R(K)~. 134#SN(K)-. 8367*SN(K)*#2 
* +6. 206*SN(K)#R(K)-2. 59 1*SN(K)*AR(K) 
CYXP(K)=5.348-6.56*R(K) ~2. 161*SN(K)-. 105#SN(K)**2 
* +45. S34HSNCK)#R(K)~. 277 *SN(K)HER(K) 
CYYP(K)=1. ~..64114R(K)~. 797 1SN(K)+. 4B45#SN(K) 12 
* +45. 3996*SN(K)#R(K)+1. 6987#SN(K)*ER(K) 
ELSE 
CXXP(K)=151. 70-220. O5#R(K)+472. 18*SN(K)-867. 94*SN(K)**2 
* #59. SO*SN(K)#R(K) -224, 268SNCK)*AR(K) 
CXYP(K)=9..272-14. 117#R(K) -16. 803*SN(K)~21. 363*SN(K)**2 
* +44, BO2*SN(K)*R(K)=1.517*SN(K)R(K) 
CYXP(K)=. 25~.0677*R(K)~64.573*SN(K)+60. 894#SN(K)##2 
* +63. 7LTASN(K)*R(K)413. 8748SN(K)*#R(K) 
CYYP(K)=10. 58-14.77#R(K)+63. 23%SN(K)-87.318SN(K)**2 
* 2B. 93*SN(K)*R(K)~18, 224SN(K)*AR(K) 
ENDIF 
ENDIF 
ENDIF 
222 CONTINUE 
NCL) 
C2) 
=XXKP(1)*WRTI/C. 


-YXKP(1)@WRTI/C 
-YXKP(2)*WRT2/C. 
YKP(1)*RTI/C 
‘YKP(2)*WRT2/C 


GXXI=CXXP(1)#NRTI/(C*W1) 
(CXX2=CXXP(2)*WRT2/(C*W2) 
CXY1=CXYP(1)*WRT1/(C#W1) 
CXY2=CXYP(2)WRT2/(CAW2) 
GYX1=CYXP(1)4WRTI/(C*W1) 
GYX2=CYXP(2)*WRT2/(C*W2) 
CYY1=CYYP(1)*WRT1/(C#W1) 
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CYY2=CYYP(2)*WRT2/(C#W2) 


RETURN 
END 


 trrrrctriiiaiiokinhdnhii fie 


ce * 
© # SUBROUTINE FOR DECISION OF THE SIGN * 
ce * 


CG etdnnnannnaiiriro irre 


SUBROUTINE SIGN(X,S) 


IMPLICIT REAL*8 (A-H) 
IMPLICIT REAL*8 (0-2) 


IF(X. GT. .0)THEN 


ELSE 
TF(X. EQ. .0)THEN 
20 


ELSE 


ENDIF 


RETURN 
END 
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aaacea 


sere sheet 


SUBROUTINE TO EVALUATE PSI'S WHICH WILL BE USED TO DEFINE 
‘THE NEXT POSITIONS ON PHASE PLANE 


SUBROUTINE PSIS(DIS, VEL, DELTA, PSI,P) 


IMPLICIT REAL*8 (4-H) 
IMPLICIT REAL*8 (0-2) 


IF(DIS#DELTA. EQ. .0)THEN 
PSI=3, 14159/2 
ELSE 
PSI=ATAN(ABS( VEL) /P/ABS(DIS#DELTA)) 
ENDIF 


RETURN 
END 


sees 


a7 


aaaaca 


tients tte 
‘ ‘THIS SUBROUTINE DEFINES THE NEXT POSITIONS * 
ON THE PHASE PLANE * 


{eeeideirkiekidkirkinkinkinkie i irre 


SUBROUTINE PHASE(DIS, VEL, DELTA, PSI, ANG, SGN1,SGN2,P, DT) 


IMPLICIT REAL*B (A-H) 
IMPLICIT REAL*8 (0-2) 


IF(DIS#DELTA. GT. .0)THEN 
IF(VEL. GT. .0)THEN 
ANG=PSI-P*DT/2. 


ANG=PSI+PADT/2. 
sis 


ELSE 
IF(DIS#DELTA. LT. .0)THEN 
IF(VEL. GE. .0)THEN 
ANG=PSI+P*DT/2. 


IF(VEL. GT. .0)THEN 
‘ANG=3.14159/2. -P#DT/2. 


IF(VEL. LT. -0)THEN 
-14159/2. -P#DT/2. 


ENDIF 
ENDIF 
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